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ABSTRACT 

The  ocean  near  surface  kinematics  were  measured  as  part  of  the  SAXON  experiment 
during  the  period  from  5  to  14  October  1988.  A  growing  wave  field  due  to  a  passing  front 
with  winds  up  to  17  ms-'  occurred  during  the  first  5  days.  A  wide  range  of  wind  forcing 
and  wave  conditions  occured  during  the  measurement  period.  Two  orthogonal  compo¬ 
nents  of  the  horizontal  velocity  were  measured  using  an  electromagnetic  current  meter. 
The  wave  directional  field  was  measured  using  a  high  resolution  2  meter  square  slope 
array.  The  classical  Longuet-Higgins  et  al  (1963)  method  of  computing  directional  wave 
spectra  and  the  new  exact  Fourier  coefficients  representation  method  (Grauzinis,  1989) 
are  used  to  compute  directional  wave  spectra.  The  new  method  of  computing  directional 
wave  spectra,  which  represents  bi-modal  distributions  of  wave  energy  exactly  matching 
the  measured  Fourier  coefficients  to  second  order,  demonstrated  improved  performance 
over  the  classical  technique.  Phase  functions  between  the  sea  surface  elevation  and  hor¬ 
izontal  velocity  components  compared  well  with  the  results  predicted  by  linear  wave 
theory.  The  coherence  function  of  the  horizontal  velocity  components  and  the  sea  sur¬ 
face  elevation  confirmed  the  dependence  on  the  wave  field  directivity.  Spectra  of  the  sea 
surface  elevation  computed  from  pressure  and  horizontal  velocity  data  compared  well 
using  linear  theory  transfer  functions.  _ _ _ 
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I.  INTRODUCTION 


A.  THE  PROBLEM 

The  need  for  the  understanding  of  the  kinematics  of  the  upper  ocean  boundary  layer 
has  long  been  recognized.  The  motions  in  the  upper  ocean  provide  the  means  for  the 
transfer  of  energy  between  the  atmosphere  and  the  deep  ocean.  The  general  ocean  cir¬ 
culation  has  its  roots  in  this  exchange,  and  the  understanding  of  the  global  circulation 
is  dependent  on  the  knowledge  of  the  upper  ocean  response  to  the  external  forcing. 
Considerable  effort  has  been  made  to  characterize  the  fundamental  mechanisms  gov¬ 
erning  the  kinematic  behavior  of  the  near  surface  ocean  boundary  layer.  However  some 
fundamental  questions,  such  as  the  nature  of  the  turbulence  below  the  ocean  surface, 
are  still  to  be  answered. 

Oceanography  has  often  followed  developments  by  meteorologists  to  improve  the 
understanding  of  ocean  processes.  However  there  are  differences  between  the  oceanic 
and  the  atmospheric  boundary  layers  that  do  not  allow  parallel  developments.  Radiative 
and  evaporative  processes,  turbulence  generated  by  breaking  surface  waves,  the  impor¬ 
tance  of  the  mean  flow,  the  dominant  signal  of  the  wave  induced  motions  and  the  dif¬ 
ferent  scale  of  the  phenomena  all  are  differentiating  factors. 

The  strong  non-linear  nature  of  the  free  surface  boundary  conditions  makes  the 
theoretical  study  of  wave  instability  and  the  onset  of  breaking  difficult.  Breaking  waves 
are  recognized  as  having  a  primary  role  in  the  transfer  of  momentum,  heat  and  gas 
across  the  air-sea  interface.  The  large  number  of  sources  of  turbulence  and  their  over¬ 
lapping  characteristic  frequency  bands  make  the  identification  of  their  relative  impor¬ 
tance  difficult.  Moreover,  the  wave  field  has  a  strong  signal  throughout  the  upper  ocean 
water  column  that  increases  the  difficulty  in  isolating  the  small  turbulence  signal.  These 
effects  place  great  demands  on  instrumentation  to  measure  the  characteristic  field  vari¬ 
ables. 

B.  THE  OBSERVATIONAL  EFFORT 

The  Naval  Postgraduate  School  (NPS)  participated  in  the  SAXON  experiment 
conducted  during  September  and  October  of  1988  off  the  East  coast  of  the  United 
States.  The  Chesapeake  Light  Tower  was  used  as  a  platform  for  investigators  and  in¬ 
struments.  The  main  goal  of  the  experiment  was  to  improve  the  understanding  of  the 
mechanisms  responsible  for  SAR  imaging  of  the  ocean  surface  and  upper  ocean  dy- 
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namics.  To  achieve  this  goal,  a  large  set  of  measurements  were  acquired  including  wind 
stress,  radar  backscatter,  directional  wave  spectra  and  near  surface  kinematics. 

The  instrumentation  deployed  by  the  NPS  group  consisted  of  a  bistatic  coherent 
accoustic  doppler  profiling  velocimeter,  a  slope  array  of  pressure  transducers,  and  an 
electromagnetic  current  meter.  The  variables  measured  included  two  orthogonal  com¬ 
ponents  of  the  horizontal  velocity  field,  directional  wave  field  and  the  three  orthogonal 
components  of  the  velocity  field  throughout  the  watei  column.  Meteorological  data  and 
wave  amplitude  information  were  measured  by  the  Massachusets  Institue  of  Technology 
(MIT).  The  bistatic  coherent  acoustic  doppler  profiler  velocity  data  are  not  reported,  as 
data  processing  methods  for  this  new  technology  are  still  being  developed. 

C.  THE  OBJECTIVES 

The  objective  of  this  thesis  is  to  measure  and  analyze  the  kinematics  of  the  near 
surface  ocean  boundary  layer,  which  is  shown  schematically  in  Figure  1.  Wind  and 
waves  are  the  components  of  the  forcing  applied  to  the  near  surface  layer.  The  kinematic 
response  is  represented  by  the  velocity  field,  that  can  be  thought  of  as  being  made  up 
of  two  components,  the  mean  and  the  fluctuating  velocities.  The  fluctuating  component 
of  the  velocity  can  be  partitioned  further  into  irrotational  and  rotational  parts.  The 
irrotational  contribution  is  classically  associated  with  the  wave  induced  orbital  motion, 
and  the  rotational  component  with  the  turbulence  field. 

In  the  second  chapter  of  this  thesis,  linear  wave  theory,  as  representative  of  the 
irrotational  component  of  the  fluctuating  velocity  component  is  reviewed,  two  methods 
of  computing  directional  wave  spectrum  are  presented,  and  the  characterization  of  the 
mean  and  fluctuating  components  of  the  velocity  field  is  considered  by  renewing  previ¬ 
ous  theoretical  and  observational  results.  The  experiment,  the  instrumentation  deployed 
and  their  operational  characteristics,  and  the  data  processing  techniques  are  described 
in  chapter  three.  Chapter  four  contains  analysis  of  the  results  for  the  period  5  to  14 
October  with  more  intensive  analysis  for  three  selected  one  hour  periods.  A  discussion 
of  the  results  is  also  included  in  chapter  four.  Finally,  conclusions  and  rccomendations 
are  presented  in  chapter  five. 
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II.  THEORETICAL  AND  OBSERVATIONAL  BACKGROUND 


A.  THE  LINEAR  WAVE  THEORY 
1.  Formulation 

Airy  or  linear  wave  theory-  has  been  successfull  in  decribing  many  observed  wave 
phenomena.  In  the  formulation  of  this  theory-,  see  for  instance  Kinsman  (1965),  the  basic 
assumptions  can  be  summarized  as 

•  Consideration  of  waves  whose  frequency  are  much  greater  than  the  earth's  rotation 
frequency  or,  co^>Q 

•  Assuming  the  fluid  as  homogen;  as  or,  co$>iX 

•  Surface  tension  effects  are  negligible  or,  co^c o 

•  Considering  the  fluid  motion  as  irrotational  or,  V  x  V  =  0  ,  which  allows  the  ve¬ 
locity  potential  to  be  defined  as  V  =  —  V0 

•  Considering  the  fluid  as  incompressible  or,  V  •  1=0 

where  co  is  the  surface  gravity  wave  frequency,  Q  is  the  earth's  frequency  of  rotation,  N 
is  the  Brunt  -  Vaisala  frequency,  u)cp  is  the  frequency  of  the  capillary  waves,  V  is  the 
water  particle  velocity  vector  and  4>  is  the  velocity  potential. 

The  combination  of  the  assumptions  of  irrotationality  and  incompressibility 
leads  to  the  basic  equation  of  linear  wave  theory,  Laplace's  equation 

V20  =  0.  [1] 


Analytical  solutions  to  this  equation  are  obtained  by  applying  the  following  boundary 
conditions: 

•  Waves  are  periodic  in  space  and  time 


•  The  bottom  boundary  is  horizontal  and  impermeable  ,  which  is  prescribed  bv 


[2] 


The  linearized  kinematic  free  surface  boundary-  condition  is  expressed  as 


cr>  I  n 

—  -  =0 


and  the  linearized  dynanuc  boundary  condition  as 


ct-TtU-*  M 

where  w  is  the  scalar  vertical  velocity,  h  is  the  local  water  depth,  tj  is  the  instanta¬ 
neous  free  surface  elevation,  g  is  the  acceleration  of  gravity  and  t  is  time. 

In  the  linearization  process,  it  is  assumed  that  the  wave  amplitude  a  is  small  when 

compared  to  the  wavelength  L  and  depth,  or  that  ~~  <^1  and  —  <^I.  These  restrictions 

L  h 

make  the  linear  wave  theory  strictly  applicable  to  small  amplitude  waves  of  infinitesimal 
steepness.  The  resulting  solution  to  the  wave  equation  is  a  single  harmonic  function. 
Invoking  the  superposition  principle,  the  instantaneous  sea  surface  elevation  can  be 
represented  as  an  infinite  sum  of  sinusoids  as  expressed  by 


cos (Kn  *x-(oni  +  e „)  [5] 

n- 1 

where  K  is  the  horizontal  vector  wavenumber,  x  =  (or, y)  is  the  position  vector,  co„  is  the 
circular  frequency  and  t„  is  the  initial  phase  angle.  The  dispersion  relationship  is  obtained 
from  the  solution  of  of  the  wave  equation  in  applying  the  surface  boundary  conditions 

oj2  =  gk  tanh(i kh).  [6] 

The  spectral  component  indice  n  is  now  dropped  for  convenience. 

The  solution  of  the  linear  wave  equation  can  be  used  to  describe  the  wave  in¬ 
duced  motions.  The  horizontal  water  particle  velocity  component  is  given  by 

f  cosh  k(h  +  z)  1  -  _ 

“fr’-r  smh kh  °<™(*  •*-«>'  +  «)  V] 


V{‘) 


CO 

-Y 


= 


n=  I 


and  the  vertical  component  given  by 


w(r)  = 


O) 


sinh  k{h  +  :) 
sinh  kh 


a  sin(A’  •  x  —  cot  +  e) 


[8] 


where  z  is  the  depth  of  observation.  Similarly  ,  the  pressure  at  a  depth  z  can  be  obtained 
from 


f  cosh  k(h  +  z)  1  -  ^ 

P{,)  ‘  cosh  kh  f  C0s('''  •  ■ 1  +  £) 


[9] 
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where  p  is  the  density.  Comparing  the  equations  [5],  [7]  and  [9],  it  is  readly  seen  that 
the  horizontal  velocity  and  the  pressure  are  in  phase  with  the  instantaneous  surface  ele¬ 
vation.  From  equation  [8]  it  is  concluded  that  the  vertical  component  of  the  velocity  is 
phase  advanced  by  90  degrees  with  the  sea  surface  elevation  and  the  horizontal  velocities 
when  considering  the  vertical  axis  oriented  positively  upwards  from  the  surface. 

Two  limiting  regions  can  be  considered,  the  deep  and  the  shallow  water  wave 
approximations.  These  simplified  solutions  correspond  to  the  consideration  of 

asymptotic  forms  for  the  hyperbolic  functions  entering  the  basic  solution  definitions. 

h 

These  extreme  cases  are  generally  identified  by  the  ratio  —  where  h  is  the  local  water 

depth  and  L  is  the  wavelength  as  computed  by  the  expressions  corresponding  to  each 

h  1 

one  of  the  cases.  For  -j->—  the  solution  is  considered  as  belonging  to  the  deep  water 


h 
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limit  and  for  “<Tq‘  s°luti°n  corresponds  to  the  shallow  water  case.  In  between 
the  full  solution  must  be  considered  and  the  solution  is  known  as  the  intermediate  water 
wave  solution. 

2.  Spectral  and  statistical  relations 

Open  ocean  data  are  invariably  collected  in  the  form  of  time  series  of  processes 
that  are  inherently  random.  Spectral  analysis  techniques  are  a  natural  first  choice  to 
study  these  data.  The  previously  presented  solution  to  the  Laplace's  equation  gives  rise 
to  simple  spectral  relationships  between  the  sea  surface  elevation  spectrum  and  velocity 
and  pressure  spectra  when  considering  the  water  column  as  a  constant  parameter  linear 
system.  Under  such  a  system,  any  two  quantities  are  connectedd  in  the  spectral  space 
via  a  transfer  function  as 


Sy(co)=  |  H(co) !  2Sx(co) 


[10] 


where  S,(«)  and  5,(co)  are  respectively  the  output  and  input  energy-density  spectra  and 
H(w)  is  the  transfer  function.  The  constant  parameter  linear  system  preserves  the  input 
frequency  at  the  output  and  modifies  the  amplitude  and  phase  of  each  component  of  the 
input  independently.  Similar  types  of  relations  can  be  deduced  for  the  cross-spectral 
quantities.  For  the  case  of  the  open  ocean  wave  field  and  considering  the  sea  surface 
elevation  spectra  as  input  and  the  velocities  or  pressure  as  output,  the  applicable  transfer 
functions  are,  for  the  horizontal  velocity  component 


/4(a>)  =  a) 


cosh  k{h  +  z) 
sinh  kh 


[11] 
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for  the  vertical  velocity  component 


Hw(a>)  =  —ico 


sinh  k(h  4-  z) 
sinh  kh 


and  for  the  pressure 


Hp{a>)  =  pg 


cosh  k(h  +  z) 
cosh  kh 


[12] 


[13] 


Spectrally  the  coherence  function  is  the  measure  of  linear  correlation  between 
two  quantities  that  can  be  thought  as  the  input-output  of  a  linear  system.  This  function 
is  obtained  from  the  measured  auto  and  cross-spectral  density  functions,  and  generally 
is  considered  in  the  form  of  squared  coherence  defined  as 


>'Lm  = 


1  Sxv(a>)  | 2 


[14] 


where  S,y{ co)  is  the  cross  spectral  density  function.  The  coherence  has  a  value  between 
0  and  1,  where  a  perfect  linear  relationship  between  input  and  output  corresponds  to  the 
value  of  1.  When  departures  from  that  limiting  value  occur,  it  can  be  hypothesised  that 
the  process  is  not  linear,  noise  is  present  in  the  measurements  or  there  are  other  inputs 
that  must  be  considered  (Bendat  and  Piersol,  1986). 

It  was  noted  in  the  formulation  of  the  linear  wave  theory  that  specific  phase 
relations  exist  between  the  sea  surface  elevation,  velocity  and  pressure.  The  relative 
phase  between  two  quantities  in  spectral  terms  is  defined  as 


cb^cj)  =  arcian 


Qxy(v) 

Cxy(oj) 


[15] 


where  <2„(oi)  is  the  quadrature  spectrum  and  CMy(co)  is  the  co-spectrum  of  the  two 
quantities. 

B.  THE  SPECIFICATION  OF  A  WAVE  FIELD 

The  adequate  description  of  the  wave  field  requires  knowledge  of  frequency  and  di¬ 
rectional  wave  spectra.  Directional  spectra  are  commonly  measured  by  pitch  and  roll 
buoy  and  multi-element  arrays,  which  represent  spectra  measured  at  a  point.  The  meth¬ 
ods  for  computation  of  the  directional  wave  spectrum  are  generally  classified  as 
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•  Model  fitting  methods  -  based  on  the  parametric  representation  of  the  spectrum 
as  the  classical  Longuet-Higgins  et  al  (1963)  method. 

•  Model  independent  methods  (Davis  and  Regier,  1977),  that  can  be  divided  into  a 
priori  and  data  adaptative  methods. 

A  priori  methods  require  in  advance  a  definition  of  criteria  to  be  verified  without  refer¬ 
ence  to  the  data  in  use.  The  data  adaptative  methods  are  based  in  a  posteriori  assump¬ 
tions  about  the  spectrum  ,  which  depart  from  characteristics  of  the  input  data.  Any  of 
these  processes  of  estimation  are  conditioned  by  the  traditional  constraints  inherent  in 
time  series  analysis  of  statistical  reliability  and  resolution. 

1.  Wave  field  description 

Sea  surface  elevation  is  generally  regarded  as  a  non-deterministic  random  wave 
field.  Statistical  techniques  must  be  used  to  define  the  descriptive  parameters  of  repre¬ 
sentative  data.  The  measurement  of  the  fundamental  properties  shared  by  pairs  of  ran¬ 
dom  records  in  the  amplitude,  time  and  frequency  domains  makes  use  of  joint 
probability,  cross-correlation  and  cross  spectral  density  functions.  Assuming  the  sea 
surface  displacement  is  a  stationary  random  function  of  space  and  time,  it  may  be  re¬ 
presented  as,  Cramer  (1962), 


>l(x,  0  =  JJ*  dA(K,  u))ei{K'x  -co,)  [16] 

where  A  =  ( A =  A  cos  0.  A)  =  A  sin  6)  represents  the  wavenumber  vector  and  6  is  the 
wave  direction.  From  the  above  definition  it  can  be  shown  that  if  (A,  to)  ¥=  (A',  a/) 

dA(K,  (jj)cTa{K',  co')  =  0  [17] 

and  that  if  (A,  o>)  =  (A',  a/) 

dA(K,  co)dA{K’,  at')  =  S(A,  (o)dKdo)  [18] 

where  S( A,  co)  represents  the  frequency  wavenumber  spectrum  and  the  tilde  denotes  the 
complex  conjugate.  The  spectral  quantity  [18]  can  be  thought  as  the  Fourier  transform 
of  the  correlation  function 


S(K,  <u)  =  — ^  f  f R(7,  r)e-^’r-Ml)d7dT 

(2ny  •>  J 


[19] 
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where  R(r ,  t)  is  the  correlation  function  defined  as 

R(r,  r)  =  E[t}(x,  l)rj(x  +7,1  +  t)]  [20] 

where  r  is  the  displacement  vector.  From  equation  [19],  it  can  be  seen  that  under  the 
assumption  of  stationarity  in  a  homogeneous  wave  field  the  frequency  wavenumber 
spectrum  is  real  and  positive  for  all  K  and  a>.  The  definition  of  the  wave  spectrum  from 
the  correlation  function,  equation  [20]  ,  is  used  to  relate  this  spectral  quantity  and  the 
variance 


J Js(R,  u)dKdu,  =  )?2(F.  t)  [21] 

where  the  variance  is  proportional  to  the  wave  energy. 

The  frequency  spectrum  can  be  directly  obtained  from  S{K ,  a>)  by  integration 
over  wavenumber 


S(o>)  =  S(K,  u)dK  =  -4-  I  R( 0,  r )eio3tdx  [22] 

J  (27r)  ^ 

and  the  wavenumber  spectrum  by  integration  over  frequency 

S(A')  =  [ S(K,  oj)dco  =  -4-  f  R{7, 0)e~lK  'r dr  .  [23] 

J  (2rr)2  J 

Using  the  dispersion  relation  as  given  by  the  linear  wave  theory,  here  considered  in  the 
deep  water  limit,  relations  can  be  deduced  between  the  wavenumber,  frequency  and 
wavenumber-frequency  spectra.  Following  Phillips  (1977),  it  can  be  deduced  from 
equation  [22]  that  the  variance  is  given  by 

,  0)  =  J  J  S(K,  6)  -^4  dude.  [24] 

The  frequency  spectrum  can  be  also  obtained  as  the  integral  of  the  wavenumber  spec¬ 
trum  over  all  directions  at  constant  wavenumber 
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S(a>)  =  ~~  f  S(K ,  d)de.  [25] 

g  J 

This  basic  set  of  relations  allows  for  the  determination  of  wave  field  properties  from  the 
observed  sea  surface  elevation  and  velocity  components. 

2.  Directional  spectrum 

The  summary7  of  the  wavenumber-frequency  characteristics  of  sea  waves  is 
contained  in  the  wave  spectrum  as  given  in  equation  [19].  This  quantity  is  simply  the 
three  dimensional  distribution  of  variance  generally  known  as  the  power  spectrum.  The 
knowledge  of  the  wave  field  at  fixed  positions  (Munk  et  al.,  1963),  permits  the  compu¬ 
tation  of  the  correlation  function.  The  definition  of  the  power  spectrum  as  given  by 
equation  [19]  leads  by  inverse  transforming  to  the  correlation  function,  that  is 

R{7 ,  t)  =  J  J  S(K,  (o)e{K  * '  ~wr)dKdco  [26] 

which  for  the  assumed  stationary  wave  field  is  an  even  function.  Consequently  the  power 
spectrum  is  real  and  symmetric  about  zero  frequency.  The  correlation  function  can  then 
be  rewritten  as  a  cosine  transform  of  the  power  spectrum 

R{ [7,  t)  =  JJs(A,  co)  cos  (A,  (D)dKd(D.  [27] 

The  correlation  function  can  now  be  redefined  as 


R(7 


S{K,  co)(  cos(cot)  cos (K  •  r )  —  sin(cot)  sin(A'  •  r  ))dKdco. 


[28] 


The  power  spectrum  can  then  be  calculated  in  terms  of  cosine  and  sine  trans¬ 
forms  of  the  correlation 


S{K,  co)  =  J [C(7,  «)  +  iQ(7,  c o)-]e~iK'rd7  [29] 

where  C(r ,  co)  and  Q(7,  co)  are  respectively  the  co-  and  quadrature  spectra  of  the  corre¬ 
lation.  Such  functions  could  be  easily  calculated  if  the  correlation  is  known  as  a  contin¬ 
uous  function  of  the  spatial  coordinate,  although  real  world  observations  lead  to  the 
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collection  of  data  in  a  small  number  of  positions.  Approximate  methods  of  computing 
directional  spectrum  from  observations  of  the  wave  field  at  discrete  positions  must  then 
be  considered. 

3.  Classical  approach 

Longuet-Higgins,  Cartwright  and  Smith  (1963)  (LCS)  following  the  suggestions 
of  Barber  (1946),  developed  a  method  of  computing  the  directional  spectrum  from  the 
motions  of  a  floating  buoy.  The  method  is  based  on  simple  realtionships  between  the 
cross  spectral  quantities  of  the  heave,  pitch  and  roll  of  the  buoy  and  the  Fourier  coeffi¬ 
cients  of  the  directional  spectrum.  The  directional  spectrum  can  be  represented  as  a 
Fourier  series 


S(J) 


*2r 


enl9S{f,  6)dd  =  an  +  ibn 


[30] 


where  6  is  the  wave  direction  and  the  coefficients  a  and  b  are  only  frequency  dependent 
and  represent  the  directional  distribution  of  energy.  Separating  the  frequency  and  direc¬ 
tional  contributions,  as  it  is  usually  done,  the  directional  spectrum  can  be  rewritten  as 

S(f,6)  =  S(f)D(f,d)  [31] 

where  D(f,  9 )  is  the  spreading  function.  D(f,  9)  corresponds  to  a  unity  area  weighting 
function  that  redistributes  the  energy  contained  in  the  power  spectrum  over  direction, 
ideally  reproducing  real  world  conditions.  The  three  measured  quantities  obtained  with 
a  buoy  can  be  directly  related  with  the  harmonic  representation  of  the  sea  surface.  The 
heave  corresponds  to  the  instantaneous  sea  surface  elevation  and  the  pitch  and  roll  are 
associated  with  the  spatial  partial  derivatives  of  the  surface  elevation.  Considering  a 
single  harmonic  of  the  sea  surface  representation 

rj(t)  =  a  cos(A'  •  x  —  an)  [32] 


the  two  horizontal  cartesian  components  of  the  sea  slope  can  be  written  as 

Vx([)  =  =  {  - iK  cos  eMl) 


[33] 


and 
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Vy{t) - ^—{-iKvnem.  [34] 

Cross  spectral  quantities  can  be  computed  for  the  instantaneous  surface  ele¬ 
vation  and  surface  slope,  and  related  with  the  Fourier  coefficients  of  the  directional 
spectrum.  These  cross  spectral  quantities  can  be  represented  as 


Considering  the 


Cm{f)  =  S{f)\lnD(f,  6)d6 

Jo 

[35] 

Q„nlf)  =  KS(j)\lr' D{f,6)  cos  6dd 

Jo 

[36] 

Qnr,(J)  =  KS(J)\2*D(f,6)  sin  6dQ 
y  Jo 

[37] 

c^{j)  =  K'stnJ7*  Dtf,  6)  co^ede 

V 

[38] 

C^{f)^K2S{f)\2rD(f,6)  sin 26d6 
y  y  4 o 

[39] 

G f)  =  A'2S(/)  [2,rZ)(/;  6)  cos  6  sin  6d6. 

*  y  Jn 

[40] 

directional  spectrum  in  terms  of  its  first  five  Fourier  coefficients,  that  is 


S[f,  6)  =  -j"  +  a\  cos  6  +  b}  sin  6  4-  a2  cos  26  +  b2  sin  26  [41] 

it  can  be  shown  that  these  coefficients  are  given  in  terms  of  the  cross  spectral  quantities, 
equations  [35  —  40],  by 

[12] 


«1 


1 

nK 


[43] 


12 


[44] 


*>\ 


1 

nK 


a2  = 


nK 2 


(Q.n  -Q.  J 


[45] 


b->  = 


nK 2 


[46] 


and  equivalently,  if  a  normalization  by  auto-spectrum  is  done,  the  Fourier  coefficients 
of  the  unit  area  spreading  function  are  determined.  This  approach  constitutes  the  clas¬ 
sical  LCS  method.  The  approximate  representation  of  the  directional  spectrum  in  terms 
of  its  Fourier  series  truncated  to  its  first  five  coefficients  leads  to  a  directional  distrib¬ 
ution  of  energy  as  given  by  a  cosine-bell  function.  Such  a  representation  is  of  limited 
resolution,  and  typically  has  large  negative  side  lobes  (Figure  2). 

To  solve  the  problem  of  the  unrealistic  negative  side  lobes,  LCS  imposed  a 
constraint  of  positivity  on  the  directional  spreading  function  which  resulted  in  a 
smoothing  effect.  However  this  approach  has  the  consequence  that  the  resolution  of  a 
distribution  narrower  than  cos4(-^-)  is  not  possible.  It  must  be  noted  that  this  artifice 
does  not  solve  the  real  problem  that  is  related  to  unresolved  harmonics  of  the  real  world 
distribution.  The  under  resolution  is  a  direct  consequence  of  the  severely  truncated 
Fourier  series  representation.  Also,  the  LCS  approach  is  not  capable  of  representing 
multi-modal  distributions  of  energy  that  are  frequently  observed  in  the  real  ocean  as  can 
be  concluded  from  the  truncated  Fourier  series  given  by  [41]  and  shown  in  Figure  2. 

The  method  just  described  can  also  be  applied  to  measurements  using  slope  ar¬ 
rays,  which  consist  of  multi-point  arrangements  of  pressure  transducers.  Typically  the 
sensors  in  the  slope  arrays  are  arranged  in  a  square.  The  real  slope  is  approximated  by 
dividing  the  differential  pressure  of  a  pair  of  sensors  by  their  separation  while  assuming 
the  sea  surface  slope  as  constant  between  each  pair  of  sensors.  The  method  of  computing 
the  directional  distribution  of  energy  then  proceeds  as  before  for  the  case  of  measure¬ 
ments  with  a  floating  buoy.  The  sea  surface  elevation  is  obtained  from  the  time  series 
of  pressure  by  applying  the  appropriate  transfer  function.  The  approximation  of  the 
surface  slope  introduces  errors  in  the  computation  given  by  (Seymour  and  Higgins.  1977) 
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DIRECTION  (DEG. ) 


Figure  2.  Typical  spreading  function  of  the  LCS  method:  Small  directional  resol¬ 
ution  and  unrealistic  negative  side  lobes. 


<5,= 


dxt 


$\n{KA/a,) 

KAp.l 

2 


ii 

dxt 


[47] 


where  a ,  =  (  cos  9,  sin  6)  Tor  i  =  1,2  and  A ,  is  the  horizontal  sensor  separation.  It  can  be 
seen  that  the  error  of  such  an  approximation  is  a  function  of  the  wavenumber  vector  and 
the  size  of  the  array,  decreasing  in  magnitude  as  these  quantities  are  reduced.  It  is  then 
natural  to  employ  arrays  as  small  as  possible,  limited  by  the  resolution  and  accuracy  of 
the  pressure  sensors. 

4.  Exact  Fourier  coefficient  representation  method 

The  LCS  approach  provides  a  wray  of  computing  the  first  five  Fourier  coeffi¬ 
cients  of  the  series  expansion  of  the  directional  spreading  function.  These  coefficients  are 
obtained  via  cross-spectral  analysis  of  triorthogonal  components  of  the  wave  field 
measurements.  The  number  of  Fourier  coefiicents  are  limited  because  of  the  measure¬ 
ment  technique  normally  used  -  floating  buoy  or  slope  array  of  four  sensors  in  a  square. 
Such  a  system  can  only  give  information  of  the  wave  field  through  second  order  terms. 
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The  coefficients  determined  can  be  represented  by  an  infinite  number  of  different  dis¬ 
tributions  whose  Fourier  expansion  is  common.  Examples  of  distributions  having  the 
same  Fourier  coefficients  can  be  seen  in  Figure  3.  It  is  then  natural  try  to  choose  a 
physically  realisable  distribution,  which  matches  the  coefficients  determined  by  obser¬ 
vation  and  more  exactly  simulates  the  true  directional  distribution. 

Grauzinis  (1989)  exploits  the  matching  of  a  set  of  Fourier  coefficients  by  differ¬ 
ent  distributions  to  develop  a  new  method  of  computing  directional  spectra.  The  previ¬ 
ously  defined  directional  spreading  function  D{6)  can  be  represented  as 


D{0)  =  { 1  +  2^0*  cos(A0fc)  +  bk  sin (A0*)]}  ~  [48] 

k 

where  ak  and  bk  are  the  unitary  Fourier  coefficients.  Such  coefficients  enter  the  definition 
of  the  polar  coefficients  c*  and  0*  as  given  by 

Ck  =  ak  +  ibk  =  ckelkd\  [49] 

Any  symmetric  unity  area  function  can  be  determined  solely  by  its  canonical  coefficients 
defined  as 


mk  =  A'(0)  cos  k6d9  [50] 

where  X{d)  is  any  arbitrary  unit  area  function  symmetric  around  0  =  0.  This  represen¬ 
tation  is  simply  the  cosine  transform  of  an  even  function.  It  is  seen  that  a  relation  can 
be  established  between  the  canonical  coefficients  and  the  Fourier  coefficients 

ak  +  ibk  =  cke‘k0"  =  mkelkp"  [5 1  ] 

where  the  polar  form  is  now  used  for  the  Fourier  coefficients  of  the  spreading  function. 
Using  linear  superposition,  any  unit  area  function  can  be  reproduced  by  a  weighted  sum 
of  symmetric  unit  area  density  functions  oriented  at  different  directions  as  given  by 


[52] 
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DIRECTION 


Figure  3.  Distributions  of  energy:  Ihe  three  distributions  have  a  common 
Fourier  expansion. 

where  wt  is  the  weighting  factor  and  /?,  is  the  orientation  of  each  of  the  components 
considered  in  the  series.  Equation  [52]  constitutes  the  generating  equation  for  the  polar 
coefficients  where  on  the  right  hand  side  wp  mt[f)  and  /?,  are  unknowns.  The  left  hand  side 
of  [52]  is  determined  from  the  cross  spectral  quantities  of  the  wave  field  typically  known 
through  second  order.  A  system  of  equations  can  be  established  as 


mx[wxelp'  +  w2e^J]  =  c]e19' 

[53] 

m2lw\e12^1  +  w2el2^]]  =  c2e129' 

[54] 

which  can  be  seen  to  be  underdetermined.  It  must  be  noted  that  in  the  system 
[53  —  54]  it  is  assumed  that  the  canonical  coefficients  are  common  to  both  of  the  dis¬ 
tributions  of  energy  considered.  T  his  was  the  approach  followed  by  Grauzinis  (1989)  to 
obtain  analytic  solutions  for  the  system  of  generating  equations.  By  the  fortunate  con¬ 
dition  of  orthogonality  of  the  exponentials  of  different  order  and  the  side  condition  that 
the  sum  of  the  weights  must  be  equal  to  unity,  analytic  solutions  are  formed.  The  ana¬ 
lytic  solutions  constitute  an  implicit  underdetermined  system  of  four  equations.  The 
unknowns  of  this  system  are  the  two  canonical  coefficients,  the  orientation  of  the  two 
unit  area  density  functions  and  one  of  the  weighting  functions.  To  solve  this  system, 
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further  constraints  are  needed.  Assuming  a  functional  relationship  between  the  canonical 
coefficients  m,  and  m2  particular  solutions  are  obtained.  The  original  relationship 
adopted  was  a  power  law  between  the  two  canonical  coefficients 

m2  =  m\  £55] 

that  for  different  values  of  r  originated  different  distributions  of  energy  as  analytic  sol¬ 
utions.  Such  distributions  are  of  decreasing  sharpness  as  the  magnitude  of  the  exponent 
increases.  For  the  values  of  the  exponent  2,  3  and  4,  the  resulting  distributions  in  its 
unimodal  form  are  represented  in  Figure  4,  for  two  different  values  of  the  canonical 
coefficient.  The  introduction  of  two  distributions  of  different  orientations  and  weights 
allows  this  method  to  represent  bimodal  distributions  of  energy.  Also  the  variable 
sharpness  of  the  model  distributions  considered  will  permit  the  matching  of  observed 
distributions  of  considerable  more  peakedness  than  the  classical  LCS  method.  For  a 
canonical  coefficient  of  .9  the  half-power  width  of  the  distribution  D4  is  12  degrees,  that 
is  about  one-third  of  the  typical  resolution  of  the  cosine-bell  distribution. 

Knowledge  of  Fourier  coefficients  of  order  higher  than  the  second  introduces 
considerable  advantage  to  this  method.  Some  of  the  constraints  can  be  dropped,  such 
as  the  functional  relationship  of  the  canonical  coefficients.  Moreover,  the  matching  of 
coefficients  of  higher  order  increases  the  resolution  of  the  method. 

The  exact  matching  of  the  Fourier  coefficients  through  second  order  of  any 
spreading  function  can  be  obtained  by  the  method  just  described.  The  problem  of  what 
criteria  to  use  for  the  selection  of  the  distribution  that  best  matches  the  real  data  must 
also  be  considered.  It  should  be  emphasized  that  there  exist  an  infinity  of  distributions 
that  can  be  considered  to  match  the  same  set  of  Fourier  coefficients.  The  problem  of 
selecting  the  best  fitting  distribution  partly  amounts  to  the  matching  of  the  ocean 
pcakedness.  Cartwright  (1963)  considered  this  problem  when  computing  directional 
wave  spectra  by  the  LCS  method.  He  concluded  that  information  could  be  obtained 
from  the  ratio  of  the  magnitude  of  the  coefficients  C*,  [49].  The  magnitude  of  such  co¬ 
efficient  is  given  by 

Q  “  \  ak  +  ■  C^6] 

Recalling  the  functional  relationship  of  the  canonical  coefficients  introduced  to  allow  for 
the  solution  of  the  generating  system,  [55]  it  is  concluded  that  in  the  Grauzinis  method 
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Figure  4.  Three  model  distributions:  Uni-modal  representation  of  three  different 
model  distributions  for  two  values  of  the  canonical  coefficient.  D2  ,  D3 
and  D4  correspond  respectively  to  the  exponent  2,  3  and  4  for  the  func¬ 
tional  relationship  of  the  canonical  coefficients  [55]  .  The  sharpness  of 
the  distribution  increases  with  the  magnitude  of  exponent.  For  a  value 
of  the  canonical  coefficient  equal  to  .9  the  half  power  width  of  D4  is  12 
degrees. 


it  is  more  natural  to  consider  the  ratio  of  the  natural  logarithms  of  the  Ck  coefficients, 
that  is 


LogRatio  = 


In  (C2) 
In(Cj) 


[57] 
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which  represents  the  greatest  power-law  of  the  canonical  coefficients  that  can  match  the 
input  Fourier  coefficients.  Cartwright  (1963)  verified  that  the  ratios  obtained  from  the 
observed  data  were  consistently  greater  than  the  values  corresponding  to  the  cosine-bell 
function  of  the  LCS  method.  It  w'as  hyphotesized  that  such  difference  could  be  due  to 
the  multimodality  of  the  wind-wave  spectrum.  Recalling  that  the  ratio  of  the  coefficients 
C*  as  an  indicator  was  deduced  from  uni-modal  considerations,  it  does  not  seem  correct 
to  make  use  of  it  to  infer  of  the  adequacy  of  multi-modal  distributions.  Grauzinis  (1989) 
points  out  that  the  natural  logarithimic  ratio  will  exceed  the  power  law  of  the  canonical 
coefficients  increasingly  as  the  beam  separation  gets  larger! .  Moreover,  the  noise  affects 
the  coefficients  Q  and  C,  differently.  Equations  [42  —  46]  show  that  the  effects  of  noise 
are  cancelled  for  the  coefficient  a2  ,  and  consequently  the  LogRatio  will  be  decreased. 

Grauzinis  (1989)  proposes  the  Beam  Separation  Index  (BSI)  as  an  indicator  of 
the  best  fitting  distribution.  This  measure  is  deduced  from  the  analytic  solutions  of  the 
generating  system  by  considering  the  canonical  coefficients  that  would  make  the  beam 
separation  zero  for  a  given  set  of  coefficients.  The  BSI  has  a  zero  value  in  the  unsepa¬ 
rated  canonically  matched  situation,  and  provides  a  relative  measure  of  canonical  fitting 
with  least  beam  separation.  From  the  definition  of  the  BSI,  it  is  concluded  that  the 
matching  of  broad  distributions  of  energy  will  result  in  bigger  values  of  the  index,  cor¬ 
responding  to  sharper  model  distributions.  The  BSI  is  a  better  indicator  of  the  best  fit¬ 
ting  distribution  than  the  LogRatio.  BSI  by  definition  contrasts  parameters  of  the  real 
data  with  model  parameters,  rather  than  exclusively  data  parameters,  as  is  the  case  for 
the  LogRatio.  However,  it  must  be  emphasized  that  the  canonical  coefficients  used  for 
the  comparison  correspond  to  constrained  analytic  solutions,  so  that  only  a  relative 
measure  is  obtained. 

C.  THE  VELOCITY  FIELD 

1.  Means  and  fluctuations 

The  ocean  is  constantly  being  fed  momentum  by  the  wind.  Although  several 
theories  exist  that  try  to  explain  the  transfer  mechanism,  the  exact  coupling  is  not 
known.  To  understand  the  way  momentum  is  injected  into  the  ocean  through  the  air-sea 
interface  and  to  look  for  its  depth  penetration  in  the  water  column,  it  is  necessary  to 
examine  and  characterize  the  upper  ocean  velocity  field.  To  gain  insight,  the  velocity 
field  is  often  partitioned  into  mean  and  fluctuating  parts  and  expressed  as 

1  The  beam  separation  is  defined  as  the  angle  between  the  orientations  of  the  two  unit  area 
model  density  functions. 
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v  =  v  +  v'  [58] 

where  v  is  the  total  observed  field,  v  is  the  mean  component  and  v'  is  the  fluctuating  part. 
The  fluctuating  part  can  be  further  partitioned  into  the  wave  and  turbulent  contributions 

v'  =  vw+v,  [59] 

where  subscripts  w  and  i  refer  to  the  wave  and  the  turbulent  components.  The  fluctuat¬ 
ing  part  is  random  in  nature  and  only  treatable  statistically.  An  objective  is  to  separate 
the  two  parts  of  the  total  field  spectrum  in  order  to  define  the  scales,  and  isolate  the 
spectrum  of  interest.  Such  separation  is  usually  done  by  an  averaging  process  that  must 
be  capable  of  maintaining  the  coherence  between  the  separated  parts  and  the  corre¬ 
sponding  equations  of  motion.  In  fact,  an  incoherent  averaging  process  can  result  in  the 
loss  of  understanding  ol  the  phenomena.  To  avoid  such  an  undesirable  situation  the 
method  must  verify  a  basic  set  of  requirements  that  can  be  summarized  (Panofsky  and 
Dutton.  1984)  as  follows 

•  The  averages  obtained  must  be  differentiable  to  any  order  as  required  by  the 
equations  of  motion 

•  The  Reynolds'  postulates  must  be  verified. 

The  Reynolds'  postulates  are  generally  presented  as 

•  The  fluctuations  must  have  a  zero  average,  or  v'  =  0 

•  The  correlation  between  the  mean  and  the  fluctuating  parts  must  vanish,  or 
v'v  =  0 

•  A  mean  quantity  must  be  unchanged  by  an  averaging  process  or,  v  =  v. 

Meeting  all  the  requirements  makes  the  averaging  process  non-trivial.  and  it  can 
be  shown  that  the  only  way  to  completely  satisfy  all  the  requirements  is  to  ensemble 
average.  Ensemble  averaging  is  not  applicable  to  the  length  restricted  series  that  consti¬ 
tute  the  generally  available  data  sets.  Other  processes  are  then  required.  Alternatives  to 
the  ensemble  average  include  time  averages  over  discrete  periods,  fitting  linear  regression 
equations  to  nonoverlapping  samples  of  the  time  series,  moving  averages  and  harmonic 
analysis  techniques.  All  these  processes  violate  to  some  extent  the  requirements  men¬ 
tioned  above,  but  as  LumJey  and  Panofsky  (1964)  have  shown,  the  situation  is  improved 
if  a  gap  is  present  in  the  spectra  of  the  phenomena. 
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2.  The  mean  velocity  field 

The  mean  component  of  the  total  velocity  field  is  at  least  in  principle  the  easiest 
to  characterize,  since  deterministic  relations  are  expected  to  completely  describe  it. 
Ekman  (1905)  pioneered  this  study.  Following  the  suggestions  of  Nansen,  a  constant 
viscosity  model  of  the  turbulent  friction  was  used.  The  solution  obtained  is  a  steady 
spiral  turn  of  the  currents  to  the  right  of  the  wind  (Northern  Hemisphere)  with  increas¬ 
ing  depth.  The  key  result  of  the  Ekman  layer  theory  can  be  presented  as 


Mr  =  ■ 


Ta  x  e, 


pf 


[60] 


in  which  the  net-transport,  Mt  ,  in  a  steady,  linear  Ekman  boundary  layer  flows  90  de¬ 
grees  to  the  right  of  the  applied  wind  stress,  r0  ,  and  has  a  magnitude  of  |  r0(pf)  1 1  re¬ 
gardless  of  the  vertical  distribution  of  stress  as  long  as  it  vanishes  at  some  depth.  In  the 
above  definition,  f  is  the  Coriolis  parameter  and  e%  the  vertical  unit  vector.  The 
parametrization  of  the  near-surface  stress  in  terms  of  an  assumed  constant  eddy  viscosity 
constitutes  the  weak  point  of  this  theory,  since  turbulence  in  nature  varies  considerably 
in  space  and  time. 

Prandtl  (1925)  introduced  the  mixing-length  theory  to  describe  turbulent  effects 
in  the  mean  flow  (see  for  instance  Holton,  1979).  The  basic  idea  is  the  parametrization 
of  the  small-scale  eddy  motions  in  terms  of  the  large  scale  mean  flow.  The  primary  hy¬ 
pothesis  is  that  in  a  boundary  layer  a  parcel  of  fluid  which  is  displaced  away  from  the 
solid  boundary  will  carry  the  mean  horizontal  momentum  of  the  original  level,  and  that 
after  moving  a  characteristic  distance,  analagous  to  the  mean  free  path  in  molecular 
viscosity,  the  parcel  of  fluid  will  transfer  the  excess  of  momentum  at  the  new  level,  cre¬ 
ating  a  turbulent  velocity  fluctuation.  The  velocity  fluctuation  is  proportional  to  the  free 
path  and  to  the  shear  of  the  mean  velocity.  In  the  boundary  layer,  the  velocity  scale  is 
given  by  the  friction  velocity  defined  by  u.  =  /  (~^~)  ■  The  relevant  length  scale  is  the 
distance  of  the  particle  from  the  solid  boundary.  The  mean  velocity  gradient  is  expressed 
as 


ft/ 

cz 


[61] 


and  when  vertically  integrated  leads  to  the  well  known  logarithmic  velocity  profile 
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[62] 


where  k  is  the  von  Karman's  constant  and  z0  is  called  the  roughness  parameter.  This  law 
holds  quite  accurately  in  constant  stress  layers. 

Bye  (1965)  made  observations  of  the  wind  driven  circulation  in  the  vertical  plane 
of  an  unstratified  lake  of  constant  depth  using  wooden  floats  of  variable  lengths.  Con¬ 
sidering  the  translational  velocity  of  the  floats  as  representative  of  the  vertically  aver¬ 
aged  current  in  which  they  drift,  measurements  were  made  in  the  first  meter  of  the  water 
column.  The  observations  showed  that  the  near  surface  currents  were  logarithmically 
distributed  in  good  agreement  with  the  conclusions  of  the  mixing  length  theory. 

Wu  (1975)  performed  measurements  of  the  wind  induced  currents  in  a  wind- 
wave  tank,  and  several  wind  conditions  were  simulated.  The  data  were  collected  by 
tracking  submerged  wooden  floats  and  also  by  using  a  Pitot-static  tube  that  was  placed 
at  selected  depths.  The  total  drift  current  was  considered  as  representing  the  wind  in¬ 
duced  current,  since  Wu  concluded  that  the  Stokes'  drift  was  negligible.  The  results  ob¬ 
tained  showed  a  logarithmic  Prandtl  type  distribution  for  the  current  near  the  surface 
while,  for  the  region  immediately  below  the  surface,  the  current  varied  linearly  with  the 
depth,  consistent  with  the  existence  of  a  viscous  sub-layer.  The  results  obtained  by  Wu 
are  presented  in  Figure  5. 

Churchill  and  Csanady  (1983)  measured  near  surface  water  currents  in  the 
coastal  zone  of  Lake  Huron  and  Cape  Cod  Bay  by  tracking  drifters  and  drogues  both 
visually  and  accoustically.  The  observed  near  surface  currents  showed  a  nearly  logarith¬ 
mic  variation  with  depth  to  about  1  meter  below  the  surface.  The  direction  of  the  cur¬ 
rents  was  almost  parallel  to  the  wind  and  predominant  wave  direction.  Figure  6,  show’s 
a  typical  profile  obtained  from  the  measurements.  For  this  experiment,  the  roughness 
parameter  was  calculated  by  relating  the  velocity  profile  projected  to  the  surface  velocity. 
The  estimated  values  of  the  roughness  parameter  W’ere  considerably  larger  than  typical 
values  on  the  air  side  for  similar  w’ind  conditions.  The  conclusions  obtained  by  Churchill 
and  Csanady  about  the  roughness  parameter  are  quite  different  than  those  of  Wu,  who 
concluded  that  the  roughness  length  was  smaller  in  the  water  than  in  the  air  in  low  wind 
velocities  and  much  smaller  under  high  wind  velocities.  For  a  similar  w’ind  condition, 
Wu  shows  values  of  the  roughness  parameter  of  the  order  of  a  few'  milimeters,  while 
Churchill  and  Csanady  computed  a  roughness  length  of  the  order  of  ten's  of  centimeters. 

I 
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Figure  5.  Mean  drift  current  profiles  after  Wu:  (a)  From  bottom  to  top  the 
currents  were  obtained  in  order  of  increasing  wind  velocity.  The  constant 
shear  is  consistent  with  the  existence  of  a  viscous  sub-layer,  (b)  Near 
surface  drift  current  relative  to  the  moving  water  surface.  Depth  scale  is 
logarithmic.  The  suggested  profile  agrees  with  a  Prandtl-Karman  type 
profile.  Adapted  from  Wu  (1975). 


Richman,  De  Szoeke  and  Davis  (1987)  used  a  string  of  vector-measuring  current 
meters  suspended  beneath  a  surface  float  to  assess  the  near  surface  shear.  Their  results 
show  that  the  observed  shear  is  consistent  with  the  results  of  the  mixing-length  theory, 


given  by  u.(kz)*'  for  a  layer  of  about  5  meters.  The  velocities  differences  between  2.5  and 
5.5  meters  depth  as  observed  and  the  values  estimated  by  the  logarithmic  profile  law  are 
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Figure  6.  Mean  velocity  profile  after  Churchill  and  Csauady:  Typical  profile  from 
drifter  and  drogue  data  of  velocity  relative  to  1.2  cm  is  plotted  against 
depth  (logarithmic  scale).  The  line  corresponds  to  a  least  squares  linear 
regression.  Adapted  from  Churchill  and  Csanady  (1983). 

shown  in  Figure  7.  A  clear  change  of  the  current  direction  with  depth  was  also  observed. 
From  a  downwind  logarithmic  profile,  the  mean  flow  evolved  to  a  rightward  turning 
profile  consistent  with  Ekman  theory  (Figure  7).  The  transition  for  the  two  regimes  oc¬ 
curs,  for  these  observations,  at  depths  that  scale  with  . 

From  the  diverse  set  of  studies  presented,  the  mean  velocity  field  in  the  upper 
ocean  mixed  layer  appears  to  be  characterized  by 


•  A  region  close  to  the  surface  where  the  velocity  variation  with  depth  is  linear,  and 
thus  consistent  with  the  existence  of  a  viscous  sublayer.  In  such  region  the  constant 
velocity  gradient  can  be  related  to  the  stress  by 


du  _ 
dz  ~  v* 


[63] 


where  v,  is  the  effective  viscosity. 

•  A  lower  region  where  the  velocity  varies  logarithmicaly  with  depth  consistent  with 
the  Prandtl-Karman  profile. 
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Figure  7.  Mean  velocities  in  the  upper  ocean  boundary  layer  after  Richman  et 
al.:  (a)  Velocity  differences  between  2.5  and  5.5  meters  depth  (solid 


line)  compared  to  the  logarithmic  profile  law  (dashed  line),  (b)  Mean 
flow  relative  to  63  meters  depth,  x  is  the  mean  wind  stress  vector.  The 
mean  flow  turns  to  the  right  of  the  wind  stress  with  depth  increase. 
Adapted  from  Richman  et  al  (1987). 


•  Below  these  upper  regions  a  layer  where  the  rightward  turning  profile  occurs  as 
predicted  by  the  Ekman  theory. 

The  scales  of  each  of  the  layers  are  strongly  dependent  on  the  wind  forcing.  The  different 
conditions  of  the  experiments  do  not  allow  for  a  clear  comparison  of  the  results. 

3.  The  fluctuating  Field 

The  fluctuating  component  of  the  total  velocity  field  is  more  complex  than  the 
mean  component.  Two  primary  reasons  lead  to  its  complexity  and  lack  of  understanding: 

•  the  great  number  of  sources  of  variability  and  the  inherent  non-linearity  of  the 
coupling  processes,  and 

•  the  difficulty  in  separating  small  turbulent  signals  in  the  presence  of  a  strong  wave 
signal 
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The  real  world  fluctuations  in  the  near  surface  can  be  associated  with  the  ocean 
surface  waves  in  the  form  of  orbital  velocities  observed  in  a  frequency  band  above 
approximately  10~‘  Hz.  Breaking  events  and  other  non-linear  wave  interactions  will  lead 
to  a  transfer  of  energy  away  from  the  principal  frequencies.  Such  a  transfer  can  be  di¬ 
rected  both  toward  high  and  low  frequencies.  Wave  interactions  result  in  associated  en¬ 
ergy  transfer  processes  that  covers  a  wide  band  of  wavenumbers.  The  Coriolis  force  can 
cause  inertial  fluctuations  of  the  currents  in  the  mixed  layer.  Internal  waves  constitute 
another  source  of  variability,  and  the  frequencies  associated  are  typically  near  10~J  Hz, 
that  is,  below  the  Brunt  —  Vaisala  frequency.  Tidal  motions  must  also  be  considered  with 
the  highest  frequency  of  about  1 2  hours  resulting  in  a  frequency  band  overlapping  with 
the  band  of  inertial  oscillations.  Besides  all  these,  one  must  account  for  the  shear  stress 
earning  eddies.  The  overlapping  of  the  frequency  bands  is  in  fact  the  big  problem  for 
the  individual  characterization  of  the  contributions. 

Particular  attention  has  been  given  to  the  kinematics  of  the  water  particles  un¬ 
der  wind  waves.  Considering  wind  waves  as  the  primary  source  of  fluctuations  reduces 
the  range  of  frequencies  to  analyze.  The  range  of  frequencies  will  be  restricted  to  the 
region  where  the  sea  and  swell  make  their  appearance,  that  is  frequencies  of  the  order 
of  about  10-'  Hz  and  above.  Classically,  a  distinction  has  been  made  between  the  fluc¬ 
tuations  associated  with  the  orbital  motions  of  the  waves  and  the  "other"  contributors, 
these  being  generally  referred  to  as  turbulent  sources. 

4.  Velocity  field  and  surface  elevation 

As  previously  noted,  linear  wave  theory  allows  the  description  of  the  wave  in¬ 
duced  motion  at  any  depth.  For  such  estimation,  the  knowledge  of  the  basic  wave  field 
characteristics  of  wave  height  H,  wave  period  T  and  the  local  water  depth  h  are  required. 
Under  the  assumptions  of  stationarity  and  uniformity,  it  is  possible  to  compute  the  auto 
and  cross  spectra  of  the  surface  elevation  and  the  velocity  components.  From  the  spec¬ 
tral  relationships,  important  conclusions  can  be  drawn. 

The  fluctuating  field  can  be  considered  as  the  superposition  of  the  wave  and 
turbulent  components.  The  wave  component  is  generally  considered  as  the  irrotational 
part  of  the  fluid  motion,  described  by  the  linear  wave  theory.  The  turbulent  component 
is  essentially  rotational  in  character.  Seitz  (1971)  concluded  that  the  total  water  particle 
motion  at  depth  could  be  considered  as  the  sum  of  the  wave  induced  motion  and  tur¬ 
bulent  motion,  with  the  turbulent  motion  being  all  that  was  not  wave  related.  This 
model  was  obtained  by  comparing  a  measured  velocity  spectrum  with  the  one  obtained 
by  directly  measuring  the  surface  wave  field. 


26 


The  spectrum  of  the  vertical  component  of  the  velocity  can  be  expressed  in 
terms  of  the  wave  and  turbulent  components  as 

SJco)  =  SwJco)  +  SWf(cu)  +  SWkW){(o)  +  SWtWw{  cd).  [64] 

Under  the  assumption  of  statistical  independence  between  the  wave  induced  and  turbu¬ 
lent  velocity  components,  their  cross  spectral  quantities  are  identically  zero  and  the 
spectrum  of  the  vertical  velocity  can  be  simplified  to 

=  S„Jco)  +  Sw((o)  [65] 

which  allows  the  separation  of  the  turbulent  component  to  be  considered  as  given  by 

5Wf(")  =  $»•(«)  -  ■$*>>)  [66] 

that  is,  defined  in  terms  of  measurable  quantities.  The  wave  induced  velocity  can  be 
obtained  from  the  measured  sea  surface  elevation  via  the  transfer  function.  Moreover, 
the  assumption  of  statistical  independence  between  the  wave  induced  velocity  and  the 
turbulent  component  allows  for  the  representation  of  the  cross  spectrum  between  surface 
elevation  and  the  vertical  component  of  velocity  as 

SwJco)  =  SwJa>)  [67] 


where  the  contribution  by  the  turbulent  component  is  now  zero.  The  coherence  function 
between  the  surface  elevation  and  the  vertical  component  of  the  wave  induced  velocity 
is  unity  under  the  assumption  of  a  constant  parameter  linear  system,  or 


[68] 


By  definition  of  the  coherence  function  for  the  total  vertical  velocity  and  sea  surface  el¬ 
evation  and  using  [65],  [67]  and  [68] 


2  5Wh(«)S^(co) 

y^(£U)  =  (SHW(  co)  +  SWr{a>))S„(<o) 


[69] 


which  can  be  simplified  to 
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[70] 


sWw(a>)  +  sW'((o ) 


Swju) 

Sw(co) 


From  [70]  it  is  noted  that  the  coherence  function  between  the  vertical  component  of  the 
velocity  and  the  sea  surface  elevation  is  independent  of  the  directional  distribution  of 
wave  energy. 

Considering  the  wave  induced  and  the  turbulent  components  of  the  velocity  as 
correlated  adds  complexity  to  the  problem.  From  the  definition  of  the  spectrum  of  the 
total  velocity  field  [68],  the  turbulent  component  is  given  by 

S„f(o>)  =  SJto)  -  SwJco)  -  SWmWt(a>)  -  SWiWm(<o)  [71] 

where  the  cross  spectral  quantities  cannot  be  directly  obtained  from  field  measurements. 
Separating  the  turbulent  and  wave  induced  parts  in  the  definition  of  the  total  velocity 
field  [75].  the  spectral  turbulent  density  function  is  given  by 

Sw,( <o)  =  Sw{co)  -  S,Jco)  -  SWwW(co)  -  SWWk(co)  [72] 

where  all  the  quantities  are  now  related  to  the  total  and  the  wave  induced  velocity  field. 
The  spectral  density  of  the  wave  induced  velocity  can  be  written  as 

$*>)«  l/4Ml2s„M  [73] 


and  the  cross  spectral  density  function  of  the  wave  induced  and  the  total  velocity  field 
as 


S*ww(u)  =  =  H„(co)Svw(co)  =  Nw((o)Swri(a>).  [74] 

The  spectrum  of  the  turbulent  velocity  component  [76]  can  now  be  rewritten  as 

SWi(co)  =  Sw(oj)  -  |  Hw(co)  1 2S„(co)  -  Hw(oj)S^(oj)  -  HJu)SWfJ(a>)  [75] 

or  equivalentely 

SWl(o>)  =  Sw{ co)  -  I  Hw (w)  1 25,(ca)  -  2Re{Hw(a>)Swn(co)}  [76] 

and  the  last  term  takes  the  form 
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Re{Hw(<o)S^{co))  =  |  Hw{a>)  I  Q^(<o). 


[77] 


Using  the  relations  now  derived,  the  coherence  function  between  the  total  velocity  field 
and  the  sea  surface  elevation  is  given  by 


_ i^(o>)  +  5^(a>)|2 _ 

+  1/41 2S,( cu)  +  2 1  Hw(a>)  I  Qwj0>))Sv(a>) 


[78] 


and  it  is  immediately  concluded  that  the  coherence  value  is  greater  than  for  uncorrelated 
wave  induced  and  turbulent  components  of  the  total  velocity  field.  For  statistical  inde¬ 
pendence  between  those  two  components,  the  cross  spectral  quantity  of  the  turbulent 
velocity  field  and  surface  elevation  will  be  identically  zero  with  all  the  other  quantities 
entering  the  definition  of  the  coherence  for  the  correlated  case  remaining  unchanged. 

The  coherence  function  between  the  horizontal  velocity  component  and  the  sea 
surface  elevation  can  be  formed  as 


2  SVkM  sv» 

y^(03)  =  SvJa>)  +  SVi(oj)  ~  SY(aj) 


[79] 


for  statistical  independence  between  the  wave  induced  and  turbulent  components.  The 
horizontal  velocity  is  defined  in  terms  of  its  cartesian  coordinates  as 

u  =  +  uWy  +  u,  [80] 


where  isotropy  is  assumed  for  the  turbulent  velocity  component.  In  [80],  w  represents 
the  wave  component.  The  directivity  of  the  waves  must  be  accounted  for  when  consid¬ 
ering  the  horizontal  velocities.  Recalling  the  definition  of  the  frequency  wave  spectrum 


S„M=  f  S„((U,  6)dQ  [81] 

and  the  transfer  function  as  given  by  the  linear  wave  theory,  the  coherence  function  for 
the  horizontal  velocity  components  and  sea  surface  elevation  can  be  written  as 
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y\ U")  = 


[82] 


I  Hu{o)  1 2  j  cos20S>J(a),  6)d6 


Hu(co)V 


cos2dSri(a),  d)dd  +  5U; 


and 


//»r 


sin20S  (co,  0)*/# 


>’L(w)  = 


I  //„(«)  l : 


sin20S„(o),  0)^0  + 


[83] 


with  /7„(<o)  given  by  [11],  It  is  evident  from  the  definitions  of  y^,(w)  and  yj  (cu)  that  they 
are  dependent  on  the  bivariate  wave  spectrum,  which  differs  from  the  coherence  function 
between  the  vertical  component  of  the  velocity  and  the  sea  surface  elevation. 

Definitions  [70].  [82]  and  [83]  show  that  a  decrease  in  coherence  will  result  as 
a  consequence  of  an  increased  ratio  of  the  turbulent  to  the  wave  component.  Such  situ¬ 
ations  will  be  particularly  true  for  deeper  positions  in  the  water  column,  since  the  wave 
induced  component  is  proportionally  reduced.  This  process  of  filtration  by  depth  is  fre¬ 
quency  selective  with  the  coherence  function  reducing  faster  in  the  high  frequency  region 
of  the  spectrum.  It  can  be  hypothesized  that  the  energy  content  of  an  observed  spectrum 
for  which  the  coherence  function  between  surface  elevation  and  velocity  shows  small 
values,  is  primarily  turbulent  when  assuming  statistical  independence  between  the  wave 
induced  and  turbulent  components  of  the  total  velocity  field,  negligible  measurement 
noise,  and  a  completely  linear  wave  field. 

Yefimov  and  Khristoforov  (1971)  have  investigated  the  behavior  of  the  coher¬ 
ence  function  of  the  velocity  field  and  surface  elevation  assuming  statistical  independ¬ 
ence  between  the  wave  induced  and  turbulent  part  of  the  total  velocity  field.  They 
considered  a  rrunimun  level  of  coherence,  and  define  that  above  the  frequency  where 
such  value  is  observed,  the  vertical  velocity  fluctuations  become  decoupled  from  the 
surface  elevation.  The  velocity  is  then  assumed  to  be  not  primarily  wave  related,  but 
turbulent  in  nature.  Observations  have  shown  that  the  transition  frequency,  defined  ac- 
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cording  to  the  above  criteria,  shifts  to  the  lower  frequency  as  depth  of  measurement 
below  the  surface  increases,  as  would  be  expected.  Such  shifting  can  be  observed  in 
Figure  8,  and  shows  the  relative  increase  of  the  turbulent  content  towards  the  shorter 
scales.  This  situation,  apart  from  the  depth  filtration  of  the  wave  induced  components, 
can  be  related  to  turbulence  convected  by  the  mean  flow,  and  as  such  is  independent  of 
locally  generated  turbulence.  When  analysing  the  horizontal  component  of  the  velocity, 
it  was  observed  that  the  coherence  values  were  consistently  smaller  than  the  values  ob¬ 
tained  when  considering  the  vertical  component.  Also,  the  transition  frequency  was 
lower  than  in  the  previous  case.  Such  smaller  values  can  be  associated  with  an  increase 
of  the  turbulent  component  relative  to  the  wave  energy,  that  has  in  this  case  a  prefer¬ 
ential  orientation.  It  can  be  hyphothesized  that  there  is  a  considerable  transport  of  hor¬ 
izontal  momemtum  by  the  vertical  component  of  the  fluctuating  field,  which  will  lead  to 
an  increase  of  the  critical  ratio  of  turbulent  wave  energy.  The  directional  spreading  of 
the  real  wave  field  must  also  be  considered  ([S2]  and  [S3]).  Yefimov  and  Khristoforov 
concluded  that  the  spectrum  of  the  vertical  velocity  could  be  considered  as  the  sum  of 
a  wave  and  a  turbulent  spectrum,  since  the  wave  induced  component  was  verified  to  be 
well  reproduced  by  estimation  according  to  linear  wave  theory. 

Thornton  and  Kraphol  (1974)  measured  the  waves  using  a  wave  staff  and  the 
two  orthogonal  particle  velocities,  vertical  and  horizontal,  using  an  electromagnetic 
current  meter.  It  was  concluded  from  the  observations  that  there  was  a  high  correlation 
between  the  surface  elevation  and  the  vertical  component  of  the  velocity  in  the  energetic 
wave  band.  The  limit  of  high  wave  energy  fluctuations  was  established  in  the  study  by 
Thornton  and  Kraphol  as  the  frequency  where  the  velocities  observed  and  theoretically 
deduced  by  linear  theory  become  significant  different.  As  observed  by  Yefimov  and 
Khristoforov  the  upper  limit  of  high  coherence  decreased  with  increasing  depth 
(Figure  9).  The  horizontal  velocity  and  the  sea  surface  elevation  showed  a  smaller  value 
for  the  coherence  ,  which  also  agrees  with  the  results  of  Yefimov  and  Khristoforov. 

The  phase  relations  predicted  by  linear  theory  between  the  vertical  and  hori¬ 
zontal  velocity  components  have  been  considered  by  Shonting  (1970).  Time  series  of 
horizontal  arid  vertical  velocity  components  of  motion  beneath  the  ocean  free  surface 
were  measured  using  orthogonally  mounted  ducted  meters.  The  cross  spectrum  of  the 
horizontal  and  vertical  components  show  (Figure  10)  disagreement  with  the  expected 
values  of  linear  wave  theory.  A  phase  difference  of  90  degrees  would  be  expected  from 
theory,  and  conscuucntly  their  co-spectrum  should  be  ideally  zero  over  the  energetic 
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Figure  8.  Coherence  function  between  the  sea  surface  elevation  and  the  vertical 
component  of  the  fluctuating  velocity  field  after  Yefimov  and 
Khristoforov:  The  coherence  decreases  in  the  high  frequency  range  and 
the  transition  frequency  co,  shifts  towards  the  low  frequency  with  depth 
increase.  Plots  refer  to  increasing  depth  from  top  to  bottom.  Adapted 
from  Yefimov  and  Khristoforov  (1971). 

band  of  frequencies.  The  values  of  the  co-spectrum  are  20  percent  of  the  magnitude  of 
the  quadrature  spectrum  indicating  a  significative  departure  from  predictions. 

Thornton  and  Kraphol  (1974)  found  measured  phase  differences  between  waves 
and  velocity  components  (Figure  11)  of  around  90  degrees  in  the  band  of  significative 
wave  energy,  which  is  consistent  with  the  linear  wave  theory.  The  phase  relations  of  the 
horizontal  velocity  and  surface  elevation  spectra  given  by  linear  wave  theory  were  simi¬ 
larly  confirmed  (Figure  1 1). 

Cavaleri  and  Zechetto  (1987)  also  focused  their  attention  on  the  phase  re¬ 
lationship  between  surface  elevation  and  the  velocity  components.  In  an  experiment 
conducted  in  the  Adriatic  Sea,  the  velocity  field  data  were  measured  using 
electromagnetic  current  meters  placed  at  right  angles  and  the  sea  surface  elevation  was 
measured  using  pressure  transducers  and  resistance  wave  gauges.  Their  results 
(Figure  12)  are  very  close  to  those  obtained  by  Shonting  (1970).  While  the  vertical  ve¬ 
locity  component  closely  followed  a  quadrature  relation  with  the  surface  elevation,  the 
horizontal  component  is  significatly  different  from  the  predicted  in  phase  condition  with 


Figure  9.  Coherence  function  between  the  sea  suface  elevation  and  the  vertical 
component  of  the  fluctuating  velocity  after  Thornton  and  Kraphol:  The 
coherence  decreases  in  the  high  frequency  range  with  depth.  Adapted 
from  Thornton  and  Kraphol  (1974). 

the  surface  elevation.  The  phase  relation  between  the  vertical  and  horizontal  velocity 
similarly  did  not  to  agree  with  the  theory. 

The  coherence  function  behavior  has  shown  very  consistent  results  in  all  the 
studies  previously  mentioned.  The  observations  point  to  a  high  coherence  between  the 
surface  elevation  and  the  vertical  velocity  component  and  a  decreased  value  when  con¬ 
sidering  the  horizontal  velocity  component.  The  frequency  limit  for  the  wave  related 
fluctuations,  estimated  from  linear  filtration,  appears  to  be  depth  dependent,  and  a  shift 
toward  the  low  frequencies  is  observed  for  increasing  depths.  The  vertical  velocity  com¬ 
ponent,  or  conversely  the  instantaneous  sea  surface  elevation,  are  well  modeled  by  the 
linear  wave  theory  in  the  high  energy  spectral  region. 


Figure  10.  Spectra  of  the  horizontal  and  vertical  fluctuating  velocity  components 


after  Shouting:  The  existence  of  a  non-zero  co-spectrum  violates  the 


quadrature  phase  relation  predicted  by  linear  wave  theory  for  the  the 


horizontal  velocity  components.  Adapted  from  Shonting  (1970). 

The  phase  studies  show  a  different  situation.  Clearly  different  results  were  ob¬ 
tained  by  the  difFerent  investigators.  Extensive  error  analysis  have  shown  no  indications 
that  could  reverse  the  conclusions  obtained.  In  all  the  studies  the  accurate  verification 
of  the  linear  theory  phase  relationships  between  the  sea  surface  elevation  and  the  vertical 
velocity  component,  can  be  noted. 


34 


Figure  11.  Phase  function  between  the  horizontal  and  the  vertical  fluctuating  ve¬ 
locity  and  the  sea  surface  elevation  after  Thornton  and  Kraphol:  Meas¬ 
ured  phase  relations  compared  well  with  linear  wave  theory  for  (a)  the 
vertical  velocity  component  of  the  fluctuating  velocity  field 
(quadrature)  and  the  sea  surface  elevation  and  (b)  the  horizontal  com¬ 
ponent  of  the  fluctuating  velocity  field  (in-phase)  and  the  sea  surface 
elevation.  Adapted  from  Thornton  and  Kraphol  (1974). 
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Figure  12.  Phase  function  between  the  vertical  and  horizontal  velocity  components 
and  the  sea  surface  elevation  after  Cavaleri  and  Zeccheto:  The  in-phase 
relation  for  the  horizontal  component  of  the  fluctuating  velocity  and 
the  sea  surface  predicted  by  linear  wave  theory  is  not  verified,  while  the 
vertical  component  of  the  fluctuating  velocity  field  is  in  quadrature  with 
the  sea  surface  elevation  as  predicted.  The  sea  surface  elevation  spec¬ 
trum  (£,)  and  the  spectrum  of  the  vertical  ( w )  and  horizontal  (w)  ve¬ 
locity  components  are  shown  in  (A).  (C)  shows  the  co  (Pn)  and  the 
quadrature  ((?„„)  spectra  of  the  vertical  and  horizontal  velocity  compo¬ 
nents.  Adapted  from  Cavaleri  and  Zeccheto  (1987). 
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III.  EXPERIMENT 


A.  LOCATION 

The  SAXON  (SAR  and  X-Band  Ocean  Nonlinearities)  experiment  was  conducted 
from  the  L'.S.  Coast  Guard  Chesapeake  Light  Tower  located  on  the  east  coast  of  the 
U.S.  14  miles  offshore  of  Cape  Henry,  Virginia  (see  Figure  13)  during  the  period  from 
14  September  to  14  October  19S8.  The  participation  of  the  Naval  Postgraduate  School 
(NPS)  was  during  the  period  from  2  to  14  October.  The  central  goal  of  the  experiment 
was  to  improve  the  understanding  of  the  mechanisms  responsible  for  SAR  imaging  of 
the  ocean  surface  and  upper  ocean  dynamics.  The  accomplishment  of  such  an  objective 
required  a  multi  disciplinary  effort,  including  oceanography,  subsurface  hydrodynamics, 
wave  hydrodynamics,  surface  boundary  layer  mechanics  and  electromagnetics.  The  data 
acquired  by  all  the  participants  in  the  experiment  included 

•  environment  conditons  -  wind,  waves  and  current  fields, 

•  wave  surface  feature  measurements  -  stereophotography,  video  photography  and 
wave  follower. 

•  surface  tension  effects  -  high  intensity  laser,  surface  film  sensor, 

•  active  microwave  equipment  -  altimeter  and  scaterometer,  stepped-frequencv 
delta- K  radar,  coherent  stepped  frequency  CW  synthetic  range  profile  Ka-band 
radar,  dual  frequency  radar,  Ka-band  and  FM-CW  scatterometer  and  vector  slope 
gauge, 

•  passive  microwave  equipment  -  4  dual,  lineraly-polarized  radiometers,  and 

•  aircraft  -  real  aperture  radar  X-band  and  svnthetic  aperture  radar  bands  L,  C  and 
X. 

Internal  waves  are  common  at  this  site  during  the  period  from  April  through  October. 
Swells  freely  propagate  from  the  open  ocean  as  the  tower  location  is  free  from  wave  re¬ 
fraction  from  nearby  islands.  The  nearby  bottom  topography  is  characterized  by  gently 
sloping  bottom. 

The  Light  Tower  (Figure  14)  is  a  four-legged,  20  meters  square  framed  structure, 
which  housed  the  experiment's  personnel  and  equipment.  The  diameter  of  each  leg  is 
about  I  meter  and  that  of  the  connecting  elements  is  0.45  meters.  The  instruments  de¬ 
scribed  here  were  installed  on  top  of  a  9  meter  tall  bottom-mounted  tower  in 
approximatelly  14  meters  of  water.  At  14  meters  depth,  linear  wave  theory  predicts  deep 
water  wave  conditions  for  waves  of  wavelength  shorter  than  28  meters  corresponding  to 
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Figure  13.  Location  of  tlie  Chesapeake  Light  Tower:  East  coast  of  the  United 
States  14  miles  ofTshore  of  Cape  Henry,  Virginia. 

periods  less  than  4  seconds.  Shallow  water  conditions  are  observed  for  waves  of  wave¬ 
length  greater  than  280  meters,  or  periods  greater  than  approximately  13  seconds.  The 
instrument  platform  was  located  approximately  40  meters  off  the  southeast  side  of  the 
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main  structure  (Figure  14)  so  that  the  measurements  are  free  of  induced  wave  disturb¬ 
ances  by  the  Light  Tower  structures  for  waves  from  the  northeast  to  the  southwest.  The 
combined  weight  of  the  bottom  mounted  tower  and  the  instrument  frame  was  over  300 
kilograms. 

B.  INSTRUMENTATION 

The  instruments  deployed  by  the  NPS,  shown  in  Figure  15,  were 

•  one,  two  component  electromagnetic  Marsh  McBirney  current  meter, 

•  four  Paroscientific  quartz  pressure  transducers,  and 

•  a  bistatic,  coherent  acoustic  doppler  profiler  which  measured  three  components  of 
the  velocity  profile  above  the  frame  (CADVP). 

The  charactheristics  of  these  instruments  are  described  below. 

1.  Electromagnetic  current  meter 

A  Marsh  McBirney  model  512  OEM  was  used  to  provide  a  reference  horizontal 
velocity  measurement  to  assist  with  the  CADVP  velocity  estimation.  The  instrument  is 
based  on  Faraday's  principle  of  magnetic  induction  in  which  a  conductor,  the  water, 
moving  in  a  magnetic  field,  generated  within  the  probe,  produces  a  voltage  that  is  pro¬ 
portional  to  its  velocity.  Two  orthogonal  components  of  the  water  velocity  field  are 
measured  in  the  range  ±  3  mr1.  The  calculated  error  is  in  the  band  +  3  percent  with  a 
resolution  determined  by  the  electronic  noise  floor.  The  measuring  unit  is  a  rubber  cov¬ 
ered  sphere  of  .04  meter  in  diameter  mounted  at  the  end  of  a  supporting  arm  of  about 
.5  meter  long  and  less  than  .01  meter  in  diameter.  The  current  meter  was  previously 
calibrated. 

2.  Pressure  transducers 

Four  Paroscientific  series  1000  Digiquartz  Intelligent  Transmitters  were 
mounted  in  a  2  meter  square  array  to  measure  directional  wave  spectra.  Each  unit  con¬ 
sisted  of  a  pressure  transducer  and  a  digital  interface  board  providing  a  digital  pressure 
output.  The  pressure  transducer  has  two  frequency  output  signals,  one  for  the  pressure 
and  the  other  for  the  internal  temperature.  The  resolution  varies  according  to  the 
measurement  interval.  For  the  choice  used,  the  manufacturer  indicates  a  pressure  resol¬ 
ution  of  10  parts  per  million  and  a  pressure  integration  time  of  0.07  seconds.  All  four 
transducers  were  triggered  simultaneously  twice  a  second  and  then  the 
pressure  frequency  counts  were  sequentially  interrogated  from  each  unit  by  a  computer. 
The  raw  pressure  counts  and  the  temperature  counts  were  converted  using  the  man¬ 
ufacturer's  scaling  equations  in  the  data  acquisition  computer.  The  manufacturer's  cali- 
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Figure  14.  The  Chesapeake  Light  Tower:  (a)  The  structure:  four-legged  ,  20  me¬ 
ters  square  frame,  (b)  Positioning  of  the  frame  of  instruments  relative 
to  the  tower  40  meters  off  the  southeast  corner  of  the  tower  in  ap¬ 
proximately  14  meters  of  water  and  at  about  5  meters  below  the  mean 
sea  surface. 

brations  were  used  and  a  comparison  with  an  aneroid  barometer  was  made  at  the  site 
of  the  experiment  while  the  frame  was  in  the  air. 

3.  Bistatic  accoustic  doppler  profiler  (CADVP) 

The  bistatic  coherent  doppler  profiler  has  a  3  degree  beamwidth  transducer 
which  emits  a  300  kHz  sinusoidal  pulse.  The  emitted  signal  ensonifies  particulates  and 
bubbles  as  it  travels  through  the  water  column.  Backscattered  energy  is  received  and 
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Figure  15.  The  frame  of  the  instruments;  A  bistatic  coherent  acoustic  doppler 
profiler  (one  transmitter  (A)  and  four  recievers  ( D )),  one  two  axis 
electromagnetic  current  meter  (B),  and  a  2  meter  square  slope  array  of 
four  quartz  pressure  transducers  (C).  Measures  are  in  meters. 

sampled  at  8  cm  range  intervals  by  the  four  radially  placed  receivers,  shown  in  schematic 
configuration  in  Figure  15.  The  doppler  shift  between  transmitted  and  received  fre¬ 
quency  is  used  to  compute  the  velocity  of  the  ensonified  volume  which  is  assumed  to  be 
moving  with  the  current  field. 

C.  DATA  ACQUISITION  SYSTEM 

The  data  acquisition  system  used  a  HP310  computer  to  sample  the  electromagnetic 
current  meter  and  pressure  transducer,  and  a  IIP330  computer  to  initialize  and  monitoi 
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the  CADVP  system.  The  HP310  computer  had  an  Infotek  16  channel,  12  bit  A  'D  con¬ 
verter  to  digitize  the  electromagnetic  current  meter  signals,  and  the  pressure  sensors  were 
linked  via  a  RS  232  communication  loop  at  19.2  Kbaud.  The  two  channels  of  the 
electromagnetic  current  meter  were  measured  with  .1  ms  time  difference  between  each 
other  and  therefore  have  negligible  time  lag.  The  time  base  for  sampling  the  current 
meter  and  pressure  sensors  was  the  computer  internal  clock.  The  data  acquired  with 
both  sets  of  instruments  was  first  recorded  in  one  hour  files  on  hard  disk  and  later 
transfered  to  cartridge  tape. 

The  CADVP  system  had  a  programable  digitizer  and  sequencer  which  sampled  the 
complex  accoustic  signals  and  sent  the  data  to  a  digital  recorder  and  to  the  HP330 
computer.  The  time  base  for  the  data  was  generated  by  a  SNAP  microcomputer. 

During  daylight  hours  a  video  camera  was  used  to  record  the  image  of  the  sea  sur¬ 
face  over  the  location  of  the  bottom-mounted  tower  for  evaluation  of  sea  conditions  and 
breaking  events.  The  time  code  used  for  identification  of  the  video  images  was  common 
to  the  CADVP. 

D.  DATA 

Data  were  acquired  during  5  to  14  October  19S8,  and  the  quantities  measured  were: 

•  two  orthogonal  components  of  horizontal  velocity  at  4.5  meters  below  the  surface 
at  a  sampling  frequency  of  8  Hz, 

•  four  water  column  pressure  readings  at  4.5  meters  depth  at  a  sampling  frequency 
of  2  Hz,  and 

•  the  three  orthogonal  components  of  the  velocity  field  through  the  water  column 
even-  S  centimeters  from  1  meter  above  the  frame  to  near  the  surface,  estimated 
at  a  4  Hz  rate. 

The  Massachusetts  Institute  of  Technology  (MIT)  group  made  available  to  all  par¬ 
ticipants  10  minute  averaged  meteorological  data.  This  data  set  includes  measurements 
of  the  wind  speed  and  direction,  air  and  sea  temperature,  and  wave  amplitude  measured 
by  an  array  of  three  wire  gauges.  The  CADVP  data  were  not  used  in  this  study,  as  al¬ 
gorithms  to  reliably  unwrap  the  coherently  estimated  doppler  shifts  are  still  under  de¬ 
velopment. 

1.  Processing 

The  velocity  and  pressure  data  were  first  checked  and  corrected  for  infrequent 
spikes  by  a  deglitching  routine.  Bad  data  sections  of  up  to  2  points  were  replaced  by 
linearly  interpolated  values.  Both  the  velocity  and  the  pressure  were  corrected  for  the 


42 


offsets  determined  in  the  pre-experiment  calibrations  and  an  assumed  atmospheric  pres¬ 
sure  of  1010  milibars  was  subtracted  out  of  the  pressure  values. 

The  determination  of  the  pressure  at  the  central  point  of  the  array  is  needed  in 
the  computation  of  directional  wave  spectra,  the  prediction  of  the  vertical  velocity  pro¬ 
files  as  given  by  linear  wave  theory,  and  for  comparison  with  the  electromagnetic  current 
meter  and  CADVP  data.  This  value  was  calculated  by  using  a  linear  interpolation  of  the 
readings  of  the  four  pressure  sensors  deployed.  The  approximation  is  of  the  form 

P(x,y)  =  -j  (P(x+I~,y)  +  P(x, y+!~))  -  ^  (P^y)  +  Pyy(x,y))  -  ...  [84] 

where  x  and  y  are  the  central  point  coordinates,  the  +  and  -  superscripts  correspond  to 
adding  and  subtracting  half  of  the  side  of  the  square  defined  by  the  four  pressure 
tranducers  to  the  central  point  coordinates,  and  Pxx  and  Pyy  are  second  order  derivatives 
in  x  and  y  respectivaly.  Considering  the  sea  surface  as  represented  by  a  monochromatic 
wave 


>l(x,y)  =  a  cos (K  •  x  —  a>t)  [85] 

leads  to  an  expression  for  the  relative  error  of  the  central  point  pressure  approximation 
of  the  form 


^  n+l^n+l  '*■' 

Rel.  Error  =Y} M 

n=l 

where  A,  is  the  side  of  the  square  array,  and  K  is  the  spatial  wavenumber.  The  relative 
error  as  a  function  of  the  frequency  is  shown  in  Figure  16,  when  considering  the  first 
10  terms  of  the  series  defined  above.  The  error  is  less  then  10  percent  for  frequencies 
smaller  than  .4  Hz,  or  waves  of  period  greater  than  2.5  seconds. 

The  depth  of  the  instrument  frame  in  the  water  column  is  required  for  the 
processing  of  the  data  and  cross  checking  of  the  results.  The  water  depth  including  the 
tidal  variation  was  obtained  by  computing  the  mean  values  of  the  pressure  records  con¬ 
sidering  an  averaging  period  of  10  minutes.  Tilting  of  the  frame  relative  to  the  horizontal 
is  of  fundamental  importance  for  future  computation  of  the  three  orthogonal  compo¬ 
nents  of  the  velocity  using  the  CADVP  data.  A  full  period  average  of  the  four  pressure 


43 


Figure  16.  The  central  point  pressure  approximation:  The  relative  error  of  the 


central  point  approximation  versus  frequency  considering  the  first  10 


terms  of  the  approximating  series  ([86])  is  less  than  10  percent  for  fre¬ 


quencies  smaller  than  .4  IIz. 


transducers  records  was  computed  to  assess  this  tilt.  The  resulting  tilting  of  the  frame 
was  at  most  .28  degrees. 

Mean  values  of  the  velocity  components  measured  with  the  electromagnetic 
current  meter  give  information  on  the  importance  of  the  convective  processes.  Averages 
of  the  two  orthogonal  horizontal  velocity  components  were  computed  for  different  av¬ 
eraging  periods. 

The  kinetic  and  potential  energy  of  the  wave  field  and  its  evolution  was  assessed 
by  computing  the  variance  of  the  horizontal  velocity  components  and  of  the  pressure 
records.  Auto  and  cross  spectra,  phase  and  coherence  functions  of  the  time  series  were 
computed  using  the  Fast  Fourier  Transform  technique  (FFT).  The  spectral  analysis  was 
based  on  one  hour  time  series  with  subsampled  FFT's  computed  over  128  second  inter- 


44 


vals.  The  resulting  ensemble  averaged  spectral  estimate  have  approximately  56  degrees 
of  freedom.  The  data  were  subjected  to  a  trend  and  mean  removal  before  being  trans¬ 
formed,  and  to  reduce  the  effects  of  spectral  leakage,  a  Hamming  window  was  used. 

The  computation  of  directional  spectra  was  done  using  both  the  LCS  and  the 
Grauzinis  methods.  The  Fourier  coefficients  of  the  wave  field  were  computed  using 
equations  [42  —  46],  The  Grauzinis  method  was  implemented  with  the  use  of  program 
TRIDEN.  Three  model  distributions  are  used  by  this  program  to  match  the  wave  field 
Fourier  coefficients,  which  include:  the  wrapped  Cauchy  distribution 


1  * 
1  —  m , 


1  —  2m,  cos  8  +  m 


?} 


[87] 


a  distribution  not  referred  to  in  the  statistical  literature  that  corresponds  to  the  analytical 
solution  of  the  generating  system  when  an  exponent  equal  to  3  is  adopted  in  the  func¬ 
tional  relationship  of  the  canonical  coefficients  ([55])  defined  by 


D3(8)  =  - 

a.  71 


1  +  2m 


cos  8  —  m? 


]  2  4 

—  2m:  cos  0  +  m , 


[88] 


and  the  wrapped  Normal  density  distribution 


D4(8)  =  — = 


1  (0+2  nn)2 


[89] 


in  o 


where  £r2  =  -21nm,.  These  model  distributions  decrease  in  sharpness  from  the 
wrapped-Cauchy  distribution  to  the  wrapped-Normal  distribution. 

The  method  was  tested  using  synthetic  data  created  by  computing  the  Fourier 
coefficients  of  known  distributions  of  energy.  The  distributions  used  were  the  box-car, 
tri-peak  (three  delta  functions  with  amplitudes  binomiallv  distributed)  and  the  twin-spike 
(two  delta  functions  of  the  same  amplitude).  The  canonical  coefficients  of  the  box-car 
function  are  defined  by 


mk  = 


sin  (kt) 
kc 


[90] 
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Figure  17.  Testing  the  Crauzinis  method  of  computing  directional  wave 
spectra:  The  model  distributions  correspond  to  equations  [87]  ,[88] 
and  [89]  .  The  test  data  (T)  is  in  (A)  a  pair  of  box-car  functions,  and 
in  (B)  a  pair  of  tri-peak  distributions.  For  the  broad  box-car  distrib¬ 
ution  the  best  results  correspond  to  the  least  sharp  model  distribution 
D4,  and  for  the  tri-peak  test  data  the  mode!  distribution  D2  gives  the 
best  results.  D4  is  not  represented  for  test  (B). 


where  c  represents  the  width  of  the  distribution  and  k  the  order  of  the  coefficient.  For 
the  tri-peak  distribution  the  canonical  coefficients  are  formed  as 


n\k  =  —  (1  -f  cos{kc)) 


[91] 


where  c  denotes  the  separation  in  direc'ion  of  the  side  delta  functions  relative  to  the 
central  peak.  The  canonical  coefficients  of  the  twin-spike  distribution  are  given  by 

mk  =  cos(Ae)  [92] 
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where  t  is  half  of  the  angular  distance  between  the  two  delta  functions  of  the  distrib¬ 
ution.  A  sample  of  the  results  of  the  tests  is  shown  in  Figure  17.  Numerical  results 
computed  for  some  synthetic  distributions  are  contained  in  Table  1.  The  computed  di¬ 
rections  for  the  energy  distribution  compare  well  with  the  orientation  of  the  input  syn¬ 
thetic  data.  It  is  noted  that  for  broad  input  distributions,  such  as  the  box-car  function, 
the  sharpest  model  distribution  gives  the  worse  results.  In  fact,  broad  input  energy  dis¬ 
tributions  are  matched  by  the  separation  of  the  uni-modal  components  of  sharp  model 
distributions,  while  for  the  narrow  input  energy  distributions  case  the  direction  of  the 
peak  of  energy  can  be  achieved  by  rotating  one  of  the  broad  uni-modal  model  compo¬ 
nents  to  the  opposite  direction.  This  last  situation  of  sharp  input  energy  distributions 
matched  by  broad  model  distributions  can  result  in  a  relocation  of  energy  in  secondary 
directions. 

The  distribution  adopted  as  representing  the  wave  field  corresponding  to  the 

data  acquired  during  the  experiment  is  a  weighted  sum  of  the  three  distributions  used 

bv  the  prosram  TRIDEN.  The  contribution  of  each  of  the  model  distributions  is 

4  3 

weighted  on  the  value  of  the  beam  separation  index.  In  the  sum,  the  weights  are  —  ,  -r- 

»  o  o 

and  —  with  the  largest  weight  corresponding  to  the  smallest  beam  separation  index 

O  ~  ^ 

value. 
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Table  1.  RESULTS  OF  TESTS  TO  THE  GRAUZINIS  METHOD  OF  COM¬ 
PUTING  DIRECTIONAL  WAVE  SPECTRUM:  Three  synthetic  dis¬ 
tributions  (DIST)  were  used  to  test  the  Grauzinis  method  of  computing 
wave  directional  spectra,  These  distributions  were  the  box-car  (BC),  tri¬ 
peak  (TP)  and  twin-spike  (TS),  for  different  partitions  of  energy  (ER)  and 
orientations  (DIR).  These  data  are  characterized  by  a  LogRatio  (LR)  near 
4.  The  three  model  distributions  (D2,  D3  and  D4)  correspond  to  equations 
[87],  [88]  and  [89],  W  represents  the  width  of  the  box-car  uni-modal  dis¬ 
tribution,  the  distance  between  the  three  delta  functions  of  the  tri-peak 
uni-modal  distribution  and  the  angle  between  the  two  delta  functions  of 
the  twin  spike  distribution.  The  Beam  Separation  Index  (BSI),  the 
canonical  coefficients  (CNC),  the  directions  (D1  and  D2)  and  the  per¬ 
centage  of  energy  (El  and  E2)  of  the  model  distributions  are  tabulated. 


DIST 

W 

ER 

DIR 

LR 

MD 

BSI 

CNC 

El 

D1 

E2 

D2 

BC 

12 

3:1 

-20  20 

4.04 

D2 

.012 

.998 

.71 

-21.6 

.29 

18.9 

D3 

.011 

.997 

.72 

-21.3 

.28 

19.0 

D4 

.010 

.993 

.75 

-20.1 

.25 

19.9 

24 

3:1 

-20  20 

4.05 

D2 

.018 

.993 

.64 

-26.0 

.36 

17.8 

D3 

.016 

.990 

.65 

-25.1 

.35 

19.9 

D4 

.011 

.974 

.73 

-21.0 

.27 

19.2 

TP 

B 

3:1 

-40  40 

3.50 

D2 

.128 

.999 

.75 

-40.0 

.25 

39.8 

D3 

.  1 2S 

.999 

.75 

-40.0 

.25 

39.9 

D4 

.128 

.999 

.75 

-40.0 

.25 

40.0 

B 

3:1 

-20/20 

4.03 

D2 

.011 

.999 

.74 

-20.3 

.26 

19.8 

D3 

.010 

.999 

.74 

-20.2 

.26 

19. S 

D4 

.010 

.999 

.75 

-20.0 

.25 

20.0 

TS 

20 

_ 

1:1 

-20/20 

4.28 

D2 

.137 

1. 

.5 

-20.0 

.5 

20.0 

D3 

.137 

1. 

.5 

-20.0 

.5 

20.0 

D4 

.137 

1. 

.5 

-20.0 

.5 

20.0 
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IV.  RESULTS 


A.  A  FULL  PERIOD  ANALYSIS 
1.  Wind  and  Waves. 

Wave  power  spectra,  wind  speed  and  direction  and  wave  amplitude  for  the  pe¬ 
riod  5  to  14  October  are  shown  in  Figure  18.  The  wind  speed  varied  from  negligible  wind 
to  a  maximun  speed  of  about  17  ms~\  The  wind  blew  generally  from  the  north.  Ex¬ 
ceptions  are  the  periods  from  the  begining  of  of  day  4  where  the  direction  was  highly 
variable  to  the  begining  of  day  7  when  a  full  clockwise  rotation  of  the  wind  direction 
occurred.  Considering  both  wind  speed  and  direction,  the  10  day  period  can  be  divided 
into  two  parts.  The  first  part  (PI)  was  during  the  time  of  increasing  wind  from  the  north 
and  the  second  part  (P2),  a  time  of  variable  winds. 

During  PI  the  wind  speed  increased,  from  4  to  17  ms-1  in  the  first  three  days, 
followed  by  a  gradual  decrease  of  the  wind  speed  to  values  of  about  1  ms _1.  During  the 
entire  period,  the  wind  blew  from  the  north,  except  in  the  final  half  day  when  the  winds 
were  light  and  had  variable  direction.  This  evolution  is  clearly  seen  in  the  representation 
of  the  wave  power  spectra  (Figure  18).  In  the  begining  of  PI,  the  peak  wave  energy  was 
observed  at  10  seconds  period.  When  the  wind  started  blowing  with  increased  speed,  the 
peak  period  of  the  waves  shifted  to  5  seconds.  A  gradual  evolution  was  then  observed 
with  increasing  wave  periods  to  about  16  seconds  when  the  wind  started  decreasing.  This 
shift  of  the  peak  of  wave  energy  is  consistent  with  the  classic  wave  model  evolution  for 
winds  of  increasing  speed.  The  rms  wave  amplitude  plot  (Figure  18)  shows  progressively 
growing  waves,  from  an  initial  amplitude  of  .5  meters  to  amplitudes  of  about  .9  meters. 
It  is  also  noted  there  was  delay  of  the  wave  field  to  the  increase  of  the  wind  speed,  and 
a  slower  rate  of  decrease  of  the  wave  amplitude  when  the  wind  began  to  slow  down. 

P2  was  highly  variable  in  wind  speed  and  direction.  An  increase  of  the  wind 
speed  from  5  to  14  m<r'  during  day  6  while  progressively  changing  direction  resulted  in 
a  weak  signature  in  the  wave  power  spectra  and  wave  amplitude  (Figure  IS).  The  in¬ 
crease  of  the  wind  speed  during  days  7  and  8  accompained  by  a  steady  direction  is  re¬ 
flected  in  wave  spectra  with  peak  of  5  seconds  period.  However  the  evolution  of  the  long 
period  wave  field  did  not  parallel  the  PI  period. 


49 


amplitude  (MIT)  are  10  minute  averages.  Day  0  corresponds  to  5  Oc¬ 


tober  at  00:00:00. 
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2.  Horizontal  velocity 

Ten  minute  averaged  time  series  of  means  and  variances  of  the  velocity  data 
measured  with  the  electromagnetic  current  meter  and  also  the  time  series  of  the  mean 
wind  speed  and  direction  are  shown  in  Figure  19.  The  u  component  of  the  mean  hori¬ 
zontal  velocity  clearly  shows  the  influence  of  the  tidal  currents.  This  result  could  be  ex¬ 
pected  from  the  current  meter  orientation,  where  the  positive  x  axis  pointed  towards  the 
entrance  of  the  Chesapeak  Bay.  The  mean  value  of  the  u  component  varied  from 
approximately  -20  to  10  cms-1.  During  P2,  a  360  degrees  rotation  of  the  wind  direction, 
and  an  increased  value  of  the  u  component  of  the  horizontal  velocity  is  noticed.  This  is 
particularly  true  when  the  wind  was  blowing  from  the  west.  The  mean  v  component  had 
a  completely  different  behavior.  Its  values  range  from  0  to  -30  cmsr1,  and  its  evolution 
closely  follows  the  wind  field.  No  tidal  signature  is  apparent.  The  variances  of  both 
components  had  a  very  similar  behavior  to  the  wave  amplitude  time  series  (Figure  18). 
It  can  be  noticed  that  there  was  considerably  greater  variance  associated  with  the  v 
component  of  the  horizontal  velocity,  particularly  during  PI.  In  this  period  the  values 
of  the  variance  for  the  u  component  range  from  30  to  180  cm2s~3  while  the  v  variance 
varied  from  20  to  300  cm3sr3. 

3.  Air  and  sea  temperature 

The  evolution  of  the  air  and  sea  temperature  from  5  to  14  October  is  shown  in 
Figure  20.  During  almost  all  the  period,  there  was  an  unstable  boundary  layer  situation 
with  the  air  colder  than  the  water.  An  exception  is  day  7.  During  all  the  period,  the  sea 
temperature  decreased  continuosly  with  values  ranging  from  18  to  21  degrees  centigrade. 
The  air  temperature  oscilated  between  9  and  21  degrees.  The  unstable  situation  associ¬ 
ated  with  the  active  wind  conditions  resulted  in  a  well  mixed  water  column.  As  a  result 
internal  waves  were  not  observed. 

B.  SELECTED  PERIOD  ANALYSIS 

Three,  one  hour  periods  were  selected  for  detailed  analysis  during  PI 

•  6  October  from  16:10  to  17:10  (PI  1), 

•  8  October  from  8:40  to  9:40  (PI 2),  and 

•  9  October  from  6:40  to  7:40  (PI 3) 

which  are  indicated  on  the  wind  speed  time  series  (Figure  18).  The  results  for  each  of 
these  one  hour  periods  include: 

•  directional  wave  spectra  computed  by  the  Grauzinis  method  and  by  the  LCS 
method  (only  the  results  for  day  6  are  shown  for  LCS), 


Figure  19.  5  to  14  October  horizontal  velocity  summary:  Means  and  variances 

of  the  two  orthogonal  components  of  the  horizontal  velocity  measured 
with  the  electromagnetic  current  meter.  The  averaging  period  is  10 
minutes.  The  wind  speed  and  direction  are  also  represented. 


•  directional  information  as  obtained  from  the  two  orthogonal  components  of  the 
horizontal  velocity  (only  the  results  for  6  October  are  shown), 

•  spectral  density  functions  of  the  two  orthogonal  horizontal  velocity  components, 
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ages  of  the  air  and  sea  temperature  show  the  atmosphere  colder  than 
the  water  which  is  characteristic  of  an  unstable  boundary  layer  situ¬ 
ation. 


•  phase  and  coherence  functions  betw  een  each  of  the  horizontal  velocity  components 
and  the  sea  surface  elevation  computed  from  the  pressure  data,  and 

•  sea  surface  elevation  spectral  density  functions  computed  from  the  pressure  and 
horizontal  velocity  measurements. 

I.  Pll 

The  first  one  hour  period  corresponds  to  the  beginning  of  the  increasing  wind 
speed.  The  direction  of  the  winds  was  variable  ,  ranging  from  the  north-northwest  to  the 
southeast.  The  mean  wind  speed  was  6  ms~'.  The  significative  wave  height  for  the  period 
was  1  meter  and  the  peak  period  of  the  wave  field  was  9.8  seconds  (frequency  .1  Hz). 
The  mean  u  component  velocity  was  3  ems'1  and  the  v  component  -12  am-',  and  the  tide 
was  at  high  water. 

The  directional  wave  spectrum  computed  by  the  exact  Fourier  coefficients  rep¬ 
resentation  method  is  shown  in  Figure  21  and  Figure  22,  and  the  directional  wave 
spectrum  computed  by  the  LCS  method  in  Figure  23.  Ihe  directional  spectrum  coin- 
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puted  by  the  first  of  the  two  methods  shows  a  range  of  directions  in  the  high  energy  re¬ 
gion  of  the  spectrum  (.06  to  .2  Hz),  from  the  east-northeast  to  the  south-southeast.  Two 
distinct  peaks  of  energy  can  be  observed.  The  first  in  the  frequency  range  from  .08  Hz 
to  .12  Hz  (swell)  with  a  mean  direction  of  east,  and  the  second  in  the  frequency  range 
.15  to  .2  Hz  (sea)  with  a  mean  direction  of  east-northeast.  An  annulus  of  energy  can  be 
observed  in  almost  all  directions  in  the  band  of  frequencies  of  significative  energy.  The 
results  obtained  by  the  LCS  method  are,  at  first  impression,  different.  Careful  observa¬ 
tion  shows  that  the  difference  is  solely  in  the  direction  pointed  by  the  center  of  the  main 
lobe.  The  LCS  results  in  a  much  smaller  directional  resolution,  typically  greater  than 
fourty  degrees,  and  the  range  of  directions  encompassed  by  the  directional  lobe  includes, 
in  this  case,  for  each  frequency  the  results  obtained  for  the  directional  spectrum  com¬ 
puted  by  the  exact  Fourier  coefficients  representation  method.  The  LCS  method  shows 
the  previously  mentioned  negative  side  lobes  of  the  model  distribution  of  energy.  Mean 
directional  information  computed  from  the  horizontal  velocity  components  (Figure  24) 
compare  well  with  the  results  of  the  exact  Fourier  coefficients  representation  method. 

The  spectral  density  functions  of  the  two  orthogonal  horizontal  velocity  com¬ 
ponents,  coherence  and  phase  functions  between  the  u  and  v  components  of  the  hori¬ 
zontal  velocity  and  the  sea  surface  elevation  and  the  spectral  density  function  of  the  the 
sea  surface  elevation  are  shown  in  Figure  25.  The  spectral  density  function  of  the  ve¬ 
locity  components  shows  an  energy  peak  at  about  .1  Hz  which  agrees  with  the  peak  of 
the  sea  surface  elevation.  The  v  component  of  the  horizontal  velocity  has  a  higher  level 
of  energy  than  the  u  component.  Recalling  the  results  of  the  directional  wave  spectrum 
and  the  electromagnetic  current  meter  orientation,  positive  v  points  toward  the  south¬ 
east,  showing  that  the  increased  level  of  energy  for  the  v  component  is  related  to  the 
wave  field  direction.  The  coherence  function  is  different  for  the  two  velocity  components. 
The  level  of  coherence  is  greater  for  the  v  component  with  values  of  .8  for  the  high  en¬ 
ergy  frequency  range.  The  coherence  values  are  about  .6  for  the  u  component  in  the 
corresponding  frequency  range.  The  phase  function  show's  both  components  of  the  hor¬ 
izontal  velocity  in-phase  with  the  sea  surface  elevation,  a  result  expected  from  the  linear 
wave  theory.  The  spectral  density  function  observed  with  the  slope  array  and  predicted 
from  the  velocity  data  compare  well.  In  the  high  energy  frequency  range  the  plot  of  the 
sea  surface  elevation  spectral  density  function  shows  a  greater  potential  energy. 

2.  P12 

The  second  one  hour  period,  8  October  from  8:40  to  9:40,  corresponds  to  the 
end  of  the  increasing  winds.  During  this  period,  the  winds  blew'  from  the  north  with  a 
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Figure  21.  Directional  nave  spectrum  for  6  October  16:00  to  17:00  computed  by  tbe 
Grauzinis  method  (3D  representation):  Two  distinct  peaks  of  energy 
can  be  observed.  The  first  in  the  frequency  range  .08  Hz  to  .12  Hz 
(swell)  and  the  second  in  the  band  of  frequency  .15  to  .2  Hz  (sea). 

mean  speed  of  15  m$-1.  The  peak  period  of  the  wave  field  was  7.5  seconds  (.13  Hz)  and 
the  significative  wave  height  was  2.3  meters.  The  mean  horizontal  velocity  (u,v)  compo¬ 
nents  were  -2  arts-'  and  -6  cuts-'.  The  tide  was  half  way  between  the  high  and  low  tide. 

The  directional  wave  spectrum  computed  by  the  Grauzinis  method  is  shown  in 
Figure  26  and  Figure  27.  The  spectrum  is  narrower  than  in  the  previous  one  hour  pe¬ 
riod,  and  the  variance  has  increased.  The  waves  associated  with  the  peak  of  energy  are 
coming  from  the  east.  As  frequency  increases ,  there  is  a  rotation  in  the  direction  of  the 
wave  field  towards  the  north.  In  the  range  of  significant  energy,  the  spectrum  spreads 
from  the  east  to  the  northeast.  The  directional  spectrum  computed  by  the  LCS  method 
(not  shown)  had  as  in  the  previous  one  hour  period  a  different  principal  c  action  for  the 
main  lobe  of  the  model  distribution,  but  as  before  the  direction  obtained  with  the 
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Figure  22.  Directional  wave  spectrum  for  6  October  16:00  to  17:00  computed  by  the 
Grauzinis  method  (top  ■view):  As  in  the  3D  representation  an  annulus 
of  energy  in  the  high  energy  frequency  band  can  be  noticed.  The  radial 
distance  corresponds  to  frequency  and  the  orientation  of  the  plot  is  as 
specified  (north  on  top). 

Grauzinis  method  was  inside  the  main  lobe.  The  directional  information  from  the 
electromagnetic  current  meter  data  (not  shown)  was  in  close  agreement  with  the  results 


Figure  23.  Directional  nave  spectrum  for  6  October  16:00  to  17:00  computed  by  the 
LCS  method:  Shadowed  areas  represent  the  negative  side  lobes  of  the 
model  distribution,  which  is  a  characteristic  of  the  method.  Increased 
negative  side  lobes  correspond  to  frequencies  containing  less  energy. 
The  directional  lobes  contain  the  range  of  directions  computed  by  the 
Grauzinis  method. 

of  the  Grauzinis  method.  During  this  period,  notice  that  the  annulus  of  energy  previ¬ 
ously  spread  over  almost  all  directions  is  reduced. 

The  summary  of  the  spectral  quantities  computed  for  the  velocity  and  sea  sur¬ 
face  elevation  is  shown  in  Figure  28.  Comparing  the  spectral  density  functions,  the  v 
component  of  the  horizontal  velocity  shows  greater  levels  of  energy  than  the  u  compo¬ 
nent.  The  v  component  spectrum  has  a  main  peak  of  energy  at  .13  Hz  (7.7  seconds)  and 
a  secondary  peak  of  energy  at  about  .18  Hz  (5.6  seconds)  that  can  be  also  observed  in 
the  sea  surface  elevation  spectral  density  function  computed  from  the  slope  array  data. 
The  coherence  function  between  the  v  component  and  the  sea  surface  elevation  shows 


57 


the  two  orthogonal  components  of  the  horizontal  velocity.  This  infor¬ 
mation  compares  well  with  the  results  computed  using  the  Grauzinis 
method. 


high  values,  .95,  in  the  range  of  frequencies  from  .8  to  .23  Hz,  that  is  in  the  range  of  high 
energy.  Lower  levels  of  coherence  between  the  sea  surface  elevation  and  the  u  compo¬ 
nent  are  observed,  with  maximun  values  of  .65  in  the  range  of  high  energy.  The  phase 
function  shows  the  u  and  v  components  of  the  horizontal  velocity  in  phase  with  the  sea 
surface  elevation.  Good  agreement  is  observed  between  the  spectral  density  function  of 
the  sea  surface  elevation  computed  from  the  pressure  data  and  the  horizontal  velocity. 
Contrary  to  the  previous  period,  there  is  a  greater  level  of  kinetic  energy.  It  is  recalled 
that  this  period  was  characterized  by  high  wind  speeds  and  increased  significant  wave 
height. 

3.  P13 

The  last  selected  one  hour  period,  9  October  from  6:40  to  7:40,  corresponds  to 
the  light  winds  region  of  PI.  T  he  observed  mean  wind  speed  was  4  ms*1  blowing  steady 
from  the  north-northeast.  The  computed  significant  wave  height  was  1.9  meters  and  the 
peak  period  of  the  wave  field  12.8  seconds  (.08  Hz).  The  mean  velocity  components  had 
values  of  4  cms~l  for  the  u  component  and  -5  cmr1  for  the  v  component.  The  tide  was 
high  water. 

The  directional  wave  spectra  computed  by  the  Grauzinis  method  is  shown  in 
Figure  29  and  Figure  30.  The  range  of  directions  is  the  narrowest  of  the  three  selected 
one  hour  periods.  The  direction  of  the  peak  frequency  is  from  the  east.  The  directional 
spectrum  computed  by  the  LCS  method  (not  shown)  had  as  before  a  different  direction 
of  the  main  lobe.  As  in  the  first  one  hour  period,  there  is  an  annulus  of  energy  spread 
over  almost  all  the  directions.  The  directional  information  computed  from  the  horizontal 
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Figure  25.  Spectral  functions  of  the  sea  surface  elevation  and  horizontal  velocity 
components  for  6  October  16:00  to  17:00:  (A)  Spectral  density  function 
of  the  sea  surface  elevation  computed  from  the  pressure  data  (solid  line) 
and  from  the  horizontal  velocity  (dotted  line).  The  two  spectra  compare 
well.  (B)  Phase  spectra  between  the  horizontal  velocity  components  and 
the  sea  surface  elevation.  An  in-phase  relation  as  predicted  by  linear 
wave  theory  is  observed  (u  solid  line  and  v  dotted  line).  (C)  Coherence 
spectra  between  the  horizontal  velocity  components  and  the  sea  surface 
elevation.  Higher  coherence  for  the  v  component  (dotted  line)  is  ob¬ 
served.  (D)  Spectral  density  of  the  two  orthogonal  components  of  the 
horizontal  velocity  (u  solid  line  and  v  dotted  line). 

components  of  the  velocity  (not  shown)  compared  well  with  the  results  of  the  directional 
spectrum  computed  by  the  Grauzinis  method. 
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Grauzinis  method  (3D  representation):  A  narrower  spectrum  than  in 
Pll  is  observed.  Only  one  peak  of  energy  at  .13hz  (7.5  seconds).  1  lie 
annulus  of  energy  is  smaller  than  in  the  one  hour  period  of  G  October 
(Figure  21). 

A  summary  of  the  spectral  quantities  of  the  velocity  components  and  sea  surface 
elevation  is  shown  in  Figure  31.  The  spectral  density  functions  of  the  horizontal  velocity 
components  show  a  greater  level  of  energy  for  live  v  component.  Contrary  to  the  previ¬ 
ous  one  hour  period  the  peak  of  energy  for  the  two  horizontal  velocity  components  is 
now  different.  The  energy  peak  of  the  v  component  agrees  with  the  peak  of  the  spectral 
density  function  of  the  sea  surface  elevation  computed  from  the  slope  array  measured 
pressures.  The  coherence  functions  between  the  horizontal  velocity  components  and  the 
sea  surface  elevation  show  a  reduced  value  for  the  u  component  of  the  velocity  with  a 
maximum  of  .47.  The  phase  relations  as  before  show  the  u  and  v  components  of  the 
horizontal  velocity  in  phase  with  the  sea  surface  elevation.  The  sea  surface  elevation 
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Figure  27.  Directional  nave  spectrum  for  8  October  8:40  to  9:40  computed  by  the 
Grauzinis  method  (top  view):  A  narrower  range  of  frequencies  and  di¬ 
rections  than  for  the  one  hour  period  of  6  October.  Radial  distance 
corresponds  to  frequency  and  the  orientation  of  the  plot  is  as  specified 
(north  on  top). 

spectral  density  function  computed  from  the  pressure  and  from  the  velocity  data  do  not 
show  as  good  agreement  as  for  the  other  two  one  hour  periods. 


Figure  28.  Spectral  functions  of  the  sea  surface  elevation  and  horizontal  velocity 
components  for  8  October  8:40  to  9:40:  (A)  Spectral  density  function 
of  the  sea  surface  elevation  computed  from  the  pressure  data  (solid  line) 
and  from  the  horizontal  velocity  (dotted  line).  The  two  spectra  compare 
well.  (B)  Phase  spectra  between  the  horizontal  velocity  components  and 
the  sea  surface  elevation.  An  in-phase  relation  as  predicted  by  linear 
wave  theory  is  observed  (u  solid  line  and  v  dotted  line).  (C)  Coherence 
spectra  between  the  horizontal  velocity  components  and  the  sea  surface 
elevation.  The  u  component  (solid  line)  shows  a  much  smaller  coher¬ 
ence  with  the  sea  surface  elevation  than  the  v  component  (dotted  line). 
(D)  Spectral  density  of  the  two  orthogonal  components  of  the  hori¬ 
zontal  velocity  (u  solid  line  and  v  dotted  line).  More  eneigy  is  observed 
for  the  v  component  of  the  horizontal  velocity. 
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Figure  29.  Directional  nave  spectrum  for  9  October  6:40  to  7:40  computed  by  the 
Grauzinis  method  (3D  representation):  A  principal  (.08  Hz)  and  a  sec¬ 
ondary  peak  of  energy  are  observed.  An  annulus  of  energy  spread  over 
almost  all  directions,  in  the  frequency  band  of  high  energy,  is  seen. 

C.  DISCUSSION 

1.  Directional  nave  spectrum 

The  computation  of  the  directional  wave  spectrum  with  the  exact  Fourier 
coeflicents  representation  method  have  shown  results  in  good  agreement  with  the 
equivalent  information  computed  by  different  methods.  The  model  distribution  chosen 
to  represent  the  bivariate  wave  spectrum,  a  weighted  sum  of  tne  three  model  distrib¬ 
utions  used  by  program  TRIDEN,  reflects  the  weakness  of  the  available  indicators  for 
the  best  fitting  distribution.  The  analysis  showed  an  annulus  of  energy  spread  over  al¬ 
most  all  directions  in  the  frequency  range  of  significant  energy.  For  very  sharp  distrib¬ 
utions  of  energy  and  if  the  uni-modal  model  distributions  are  too  broad,  the  matching 
of  the  input  Fourier  coefficients  is  done  by  moving  one  of  the  components  away  from 
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Figure  30.  Directional  wave  spectrum  for  9  October  6:40  to  7:40  computed  by  the 
Grauzinis  method  (top  view):  The  energy  spread  over  the  west  half 
plane  has  the  largest  value  of  the  three  one  hour  periods  considered  for 
short  term  analysis.  Radial  distance  represents  frequency  and  the  plot 
is  oriented  north  on  top. 

the  component  centered  on  the  direction  of  the  energy.  The  fraction  of  energy  of  the 
relocated  uni-modal  component  will  be  minimal  but  of  increasing  value  as  the  sharpness 
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components  for  9  October  6:40  to  7:40:  (A)  Spectral  density  function 
of  the  sea  surface  elevation  computed  from  the  pressure  data  (solid  line) 
and  from  the  horizontal  velocity  (dotted  line)  do  not  compare  as  well 
as  in  the  one  hour  periods  of  6  and  8  October.  (B)  The  phase  spectra 
between  the  horizontal  velocity  components  and  the  sea  surface  ele¬ 
vation  verifies  the  in-phase  relation  predicted  by  the  linear  wave  theory. 
(C)  Coherence  spectra  between  the  horizontal  velocity  components  and 
the  sea  surface  elevation,  shows  a  significant  smaller  value  for  the  u 
component  (solid  line).  (D)  The  spectral  density  of  the  two  orthogonal 
components  of  the  horizontal  velocity  (u  solid  line  and  v  dotted  line) 
shows  different  peaks  for  the  two  spectra,  and  more  energy  is  associated 
with  the  v  component. 
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of  the  model  distribution  decreases.  This  process  is  illustrated  in  Figure  32  for  the  period 
of  8  of  October  from  8:40  to  9:40.  The  directional  spectra  shown  correspond  to  the 
wrapped-Cauchy  (D2)  and  wrapped-normal  (D3)  model  distributions,  the  first  being 
considerably  sharper.  The  amplitude  of  the  annulus  of  energy  for  D4  is  twice  the  am¬ 
plitude  for  D2.  This  relocation  of  energy  amounts  to  a  decrease  of  the  spectral  value  in 
the  range  of  frequencies  of  significant  energy.  A  reduction  of  the  error  introduced  by  this 
process  can  be  accomplished  by  considering  sharper  model  distributions.  It  is 
hyphothesized  that  the  energy  contained  in  this  annulus  is  an  indicator  of  the  best  fitting 
model  distribution.  It  should  be  noted  that  for  the  selected  one  hour  periods,  the  beam 
separation  index  for  the  majority  of  the  frequencies  was  smaller  for  the  wrapped-Cauchy 
model  distribution.  The  analysis  of  the  results  does  not  show  a  second  annulus  of  energy 
that  can  be  observed  at  the  maximun  frequency.  This  outer  annulus  of  energy  is  not  a 
consequence  of  the  method  used  to  compute  directional  spectrum.  It  represents  the  ef¬ 
fects  of  noise,  which  is  related  to  the  approximation  of  the  instantaneous  pressure  at  the 
central  point  of  the  slope  array,  and  to  the  computation  of  the  sea  surface  elevation  with 
the  use  of  the  transfer  function  given  by  linear  wave  theory,  which  amplifies  noise 
strongly  with  increasing  frequency.  The  significance  of  the  results  using  the  exact  Fourier 
representation  method  must  be  interpreted  recalling  that  only  the  Fourier  coefficients 
through  second  order  are  matched,  that  at  most  bi-modal  distributions  of  energy  can  be 
resolved  and  that  there  is  a  significant  dependence  of  the  results  on  the  relative  sharpness 
of  the  input  and  model  distributions. 

2.  Coherence  function 

The  results  obtained  for  the  coherence  functions  between  the  sea  surface  ele¬ 
vation  and  the  two  orthogonal  components  of  the  horizontal  velocity  verify  the  theore¬ 
tical  results.  In  the  theoretical  background,  it  was  noted  that  the  coherence  function 
between  the  sea  surface  elevation  and  the  horizontal  components  of  the  velocity  de¬ 
pended  on  the  bivariate  wave  spectrum.  The  results  of  the  experiment  show  that  the 
values  of  the  coherence  function  are  greater  for  the  component  that  had  the  major  con¬ 
tribution  of  the  wave  field,  that  is,  the  coherence  function  had  greater  values  for  the 
component  that  had  the  smaller  turbulent  to  wave  induced  ratio  assuming  isotropy  of 
the  turbulent  field.  An  attempt  was  made  to  compute  the  coherence  function  using  the 
information  given  by  the  directional  wave  spectrum  and  the  theoretical  definition  of  the 
coherence  function  for  statistically  independent  wave  induced  and  turbulent  components 
of  the  velocity,  equations  [82]  and  [83]  .  The  estimate  of  the  turbulence  spectrum  was 
obtained  by  linear  filtration  assuming  isotropy  of  the  turbulent  field. 
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Figure  32.  Directional  spectrum  computed  by  the  Grauzinis  method:  The  annulus 


of  energy  in  the  representative  frequency  band  is  considerably  reduced 


for  the  sharper  model  distribution  D2. 


3.  Phase  spectrum 

The  computed  phase  spectra  between  the  horizontal  velocity  components  and 
the  sea  surface  elevation  were  similar  to  the  observations  of  Thornton  and  Kraphol 
(1974)  and  the  results  of  linear  wave  theory.  For  the  three  one  hour  periods,  the  com- 
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ponents  of  the  horizontal  velocity  and  the  sea  surface  elevation  were  in-phase,  within  the 
predicted  phase  shift  due  to  the  non-centered  position  of  the  electromagnetic  current 
meter. 

4.  The  spectrum  of  the  sea  surface  elevation 

The  comparison  of  the  spectra  of  the  sea  surface  elevation  from  the  measure¬ 
ments  of  pressure  and  from  the  horizontal  velocities  show  good  agreement  for  the  three 
one  hour  periods  analyzed.  The  results  computed  for  9  of  October  from  the  two  data  sets 
are  not  statistically  different  at  the  95%  confidence  level.  It  must  be  noticed  that  the 
comparison  of  the  two  spectral  density  functions  departs  from  the  linear  wave  theory  for 
both  spectra.  The  direct  observation  of  the  sea  surface  elevation  would  introduce  more 
meaning  to  the  comparison  and  allow  for  a  significant  verification  of  the  linear  wave 
theory.  The  same  can  be  said  if  the  vertical  component  of  the  velocity  field  was  available. 
The  vertical  component  is  not  affected  by  directional  effects  as  happens  with  the  hori¬ 
zontal  velocity  components,  and  this  effect  has  been  shown  to  be  related  both  to  the 
wave  and  mean  flow  field. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  new  method  of  computing  directional  wave  spectrum  developed  by  Grauzinis 
(1989)  verifies  a  necessary,  but  not  sufficient,  condition  for  the  estimation  of  high  resol¬ 
ution  bivariate  wave  spectrum  since  only  the  Fourier  coefficients  through  second  order 
are  matched.  The  ambiguity  of  the  detailed  shape  of  the  directional  wave  spectrum  is  due 
to  the  lack  of  a  strong  indicator  for  the  best  model  fitting  distribution.  The  results  ob¬ 
tained  with  this  new  method  have  given  good  indications  of  its  capabilities.  The  possi¬ 
bility  of  resolving  bi-modal  energy  distributions  without  a  priori  or  a  posteriori 
assumptions  constitutes  a  significant  improvement. 

The  velocity  field  measurements  confirmed  previous  observations  and  theoretical 
results.  Phase  relations  show  a  behavior  as  expected  from  linear  wave  theory.  Coherence 
functions  emphasized  the  directional  effects  of  the  wave  field.  Statistically  equivalent 
results  were  obtained  for  the  spectral  density  functions  of  the  sea  surface  elevation 
computed  from  the  pressure  and  horizontal  velocity  measurements  using  the  transfer 
functions  derived  from  the  linear  wave  theory. 

The  instrumentation  deployed  during  the  SAXON  experiment  provided  a  fairly 
complete  set  of  measurements.  Improvements  can  be  considered.  The  direct  measure¬ 
ment  of  the  sea  surface  elevation  is  of  great  importance,  particularly  at  the  central  point 
of  the  frame  of  instruments.  This  measurement  gives  the  possibility  of  computing  higher 
order  Fourier  coefficients  of  the  sea  surface  and  improving  the  directional  spectrum 
computation.  When  bistatic  acoustic  doppler  profiler  data  are  available,  more  significant 
tests  of  the  spectral  relations  between  measured  quantities  will  be  possible,  and  a  far 
more  detailed  view  of  stress  transfer  and  the  turbulent  field  should  be  available. 

The  analysis  of  the  experiment  has  shown  a  wide  range  of  forcing  conditions  that 
will  allow  different  kinematic  responses  of  the  upper  ocean  boundary  layer  to  be  studied. 
The  availability  of  the  three  components  of  the  velocity  field  throughout  the  water  col¬ 
umn  will  give  a  new  meaning  to  the  observations.  A  cross  check  of  measurements  and 
methods  will  be  then  possible,  and  a  significative  improvement  on  the  understanding  of 
the  upper  ocean  boundary  layer  is  expected. 

The  review  of  theoretical  results  has  shown  the  necessity  for  new  developments.  The 
separation  of  the  strong  wave  signal  from  the  small  background  turbulence  will  provide 
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a  continuing  challenge.  The  study  of  the  directional  spreading  effects  of  the  wave  field 
in  the  signal  at  depth  also  needs  to  be  considered. 
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